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Journal of Research of the National Bureau of Standards 


Vol. 55, No. 6, December 195: Research Paper 2633 


The First Spectrum of Barium, Bal 


Henry Norris Russell ' and Charlotte E. Moore 


The analysis of the Bar spectrum is presented. 


Although a monograph based on 


homogeneous observations is needed, this paper has been prepared from existing references, 
to present unpublished analysis by the senior author and, also, because Bar is important in 
the history of atomie spectra. The first regularities attributed to two excited electrons were 


found in Bat as well as in Catand Sr. 

Three tables are included: 
observed and all classified lines; 
combinations 


A complete term list; a complete line list containing all 
a table giving the estimated intensities of the observed 
There are approximately 390 classified lines 


The limit, derived from well-established series, is 42032.4 kK, giving an ionization poten- 


tial of 5.210 volts 


The Bat spectrum occupies an important positioA 
in the history of atomic spectra. Although regular 
series terms were reported by Saunders, Paschen- 
Gétze, and Fowler (1920 to 1922) [1] the first regu- 
larities attributed to two excited electrons were 
found in Bat, as well asin Cat and Sri, in 1925 [2]. 

Shenstone and Russell suggested some revisions in 
the early analysis and discussed the perturbed series 
of Bat in 1932 [3]. They adopted as the ionization 
limit, i. e., the absolute value of the ground state, the 
value 42032.4 K, giving an ionization potential 
5.210 volts. 

Although the spectrum has been studied thor- 
oughly, a monograph based on a complete homogene- 
ous line list is needed. In view of the fact that this 
does not appear to be forthcoming, a complete line 
list and multiplet array have been prepared from the 
existing literature references, and from unpublished 
analysis by the senior author. Uniform notation has 
been adopted, and the decimals of all energy levels 
have been revised. Sullivan and Burns (table 2, 
ref. C) improved the values of a number of energy 
levels from their interferometric measurements of 
selected lines. Starting with these three-place 
values, the rest were revised from a study of all 
observed combinations. 

The complete term list is given in table 1, 
is self-explanatory. The terms are listed 
following order: 


which 
in the 


(1) Terms due to the equivalent electron pairs 

is*, 6p*?, and 5d? so far as known. 
(2) The regular series of triplet and singlet terms 
having as limit the ground term in Ba 116s *Soi;, given 
in the order of the running electrons ns, np, nd, and 
nf. 

(3) The triplet and singlet terms produced by 
two excited electrons, having as limit the term 5d ?D 
in Bau, arranged in the same order of running 
electrons. The designations for these terms have 
primes attached to the running electrons to denote 
the higher limit term in Ba tm. 


Princeton University Observatory, Princeton, N. J 
2 Figures in brackets indicate the literature references at_thefend of this text. 
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B. E. Moore observed the Zeeman effect in Bat 
in 1908 [4] and, so far as is known to the present 
authors, his work has never been superseded. His 
measurements of AX/M have been divided by the 
correction factor 1.09, a factor adopted to reduce 
them to the theoretical Landé patterns for a very 
few well-known lines. For unresolved patterns the 
formulas for blends given by Shenstone and Blair [5] 
have been used to determine observed g-values. 
The results obtained from Moore’s observations are 
as follows: 

Desig Obs. ¢)| Desig Desig 


Obs. ¢ Desig Obs. ¢ 


op 5d 'D < 4/ 6p’ 
Hp 
op 


tip; 2 


op 
7p 


Table 2 contains all of the observed and classified 
lines of Bar. The letter in column 1 indicates the 
literature reference from which the wavelength in 
column 2 is quoted. The bibliography of these 
references is given at the end of this table. Two sets 
of intensities are quoted in column 3. The left-hand 
entry 1s from the reference source indicated in 
column 1, except for those in parentheses, which are 
arc intensities estimated by Exner and Haschek 
(ref. G). The intensities in parentheses are entered 
for the lines observed with the interferometer by 
Burns and Sullivan. The right-hand entries in 
column 3 are the are intensity and temperature class 
estimated by A. S. King [6]. Brackets denote King’s 
estimated intensity in the furnace spectrum, in cases 
where this spectrum gives better resolution than his 
are spectrum. Column 4 contains the vacuum wave 
number in em™'!, labeled K for kayser. It is followed 
by AK, the difference between the observed wave 
number and the wave number calculated from the 
energy levels. The last column gives the multiplet 
designation in the notation of table 1. There are 
approximately 390 classified lines in table 2. 
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The intensities of the observed combinations in 
Bat are recorded in table 3. This table is arranted 
in two sections: (1) the multiplets produced by tran- 
sitions from the low even terms 6s? 'S, 5d? °P, 5d? 'D; 
5d ‘DD, and Sd ‘[) more precisely 6s 5d 3 1])). to 
higher odd terms, the latter being given in the same 
order as in table 1; (2) multiplets representing 
transitions from 6p*'P° and 6p’ *!P°, 6p’ *'D 
6p’ *'F° to higher even terms as given in table 1. 
All observed combinations are reported in table 3 
except the following four, which are omitted to save 


space 


Desig Ir 
t Dy-4f FG Zt 
fd (D4 1 »m 
S-7p  °P Hi) 
S-7pl l 


The intensities in table 3 are from column 3 of table 2. 
King’s estimated are intensity (or furnace intensity 
in brackets) and temperature class are entered when 
known. All intensities not by King are given in 
parentheses. For some lines no estimated intensity 
is available. In such * has 
adopted to denote that the combination has been 


cases, a ~~ ) been 
observed 

The present summary of our knowledge of the Bat 
spectrum does not serve to replace the monograph 


TABLE 2 Observed 


, 
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needed for this spectrum. Further observations of 
the Zeeman effect are desirable, “With modern 
infrared detectors the long-wave region beyond the 
photographic range offers an attractive and worth- 
while program. In particular, the even term 5d ?°F 
remains to be found. The short-wave region like- 
wise needs to be reobserved, and many of the early 
measurements throughout the spectrum could be 
improved 


The authors are particularly indebted to W. F. 
Meggers for the use of his precise observational data 
in the photographic long-wave region, in advance of 
publication. 
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Compressibilities of Long-Chain Normal Hydrocarbons 
C. E. Weir and J. D. Hoffman 


Measurements were made of the compressions of n-CisHae, n-CaHa, 
in the erystalline state at 21° C 
Within the experimental error, which was fairly large because 


n-( sH . and n-C H, 
1,000 to 10,000 atmospheres, 


n-CogH so, n-CogH 5,, 
The pressure range used was from 


of the small size of sample, no significant variation of the compression with chain length 


was observed 
by the equation 


AV/V 


0.1717 log 


The average compression of all the hydrocarbons studied may be represented 


0.6272 


2,500+ P 


he average compressibility of the hydrocarbons extrapolated to 1 atmosphere was calculated 


to be 29.8 10-° atm > + bee © 


1. Introduction 

Considerable interest has been exhibited in the 
properties of the normal long-chain hydrocarbons, 
as indicated by numerous reports on studies of these 
materials. As pointed out earlier [1],' there is need 
for data on the compressibilities of these materials, 
the only data available on the erystalline materials 
appearing to be those of Miller [2] obtained from 
X-ray measurements on n-C.,Hy and n-C.,H, 

Data on the liquids containing between 5 and 10 
carbon atoms were reported by Bridgman [3], and 
Cutler, and Webb [4] have recently 
reported measurements on n-CyH.,. n-C),;H a, and 
n-C Hyg in the liquid state. 

In the present report, data are given for com- 
pressions at 21° © for some of the ervstalline normal! 
hydrocarbons containing between 18 and 30 carbon 
atoms. 

Because of the small amount of each hydrocarbon 
available, the errors of measurement were necessarily 
rather large. Accordingly, the present results are 
of a somewhat exploratory nature, and will probably 
be subject to slight revision when data on much 
larger sumples are obtained. 


Schiessler, 


2. Experimental Method 


It was considered desirable to measure the com- 
pressibilities at low pressures, i. e., 1 to 1,000 atm, 
but efforts to do this with the present samples in 
this pressure range, using piezometric techniques, 
were The errors in such measure- 
ments were so large as to render the data of dubious 
value These errors were probably the result of 
interfacial effects arising from the fact that the only 
confining liquids considered inert toward the speci- 
mens were water and mercury— both nonwetting 
liquids. After vacuum impregnation of the speci- 
mens with either liquid, unfilled capillary 
were visible. Presumably these were filled at ele- 
vated thus giving rise to the erratic 
compressions observed. Attempts to fill the voids 
permanently at 10,000 atm were unsuccessful be- 
cause the visible capillaries emptied partially in an 
irreproducible manner on reduction of the pressure. 


I ire n bracket nel te the terature references at the end of tl paper 


unsuccessful 


spaces 


pressures, 


It was decided, therefore, to measure the ecom- 
pressions at high pressures, where errors arising 
from void spaces are negligible. The basic tech- 


niques and equipment used in such studies have 
been described earlier [5]. 
Certain modifications of the 
devised to permit study of the small specimens. 
Each hydrocarbon was fused and recrystallized in a 
small Pyrex tube, open at one end. The filled tubes 
were placed, open end down, in a. stainless-steel 
container, open at one end. <A steel spring under 
compression served to hold the open end of the glass 
with the bottom of the steel con- 


procedures were 


tube in contact 


tainer. This assembly was evacuated for ' hr, 
using a mechanical pump, and the stainless-steel 
container filled with clean mercury while under 


vacuum. To minimize danger of contamination of 
the pressure vessel with mercury, the mercury-filled 
assembly was placed in a somewhat larger telescoping 
steel container designed to catch any mercury 
inadvertently spilled from the inner container. This 
unit, consisting of telescoping tubes, glass, mercury, 
spring, and hydrocarbon, represented the compres- 
sion specimen and was immersed in a light petroleum 
distillate (Varsol) in the pressure For a 
comparison standard the same assembly was used, 
except that the glass tube contained no hydrocarbon. 
The comparison material for the hydrocarbon, there- 
fore, consisted of a similar volume of mercury. 

Compression data were calculated by means of a 
generalization of the equation of Adams [6], as 
follows: 


vessel. 


S 
k—k’ (aL | a P)+-Al’a(P—P" 
n=m-1 7 _\" Vk—k 
ka oP) 2 s |i yt 
Vast 
> k,)( - ) 


P denotes pressure in atmospheres and V, volume; 
subscript zero denotes 1 atm, and subscripts 7 and 
refer to measurements made with reference ma- 
terial and pressure-transmitting liquid, respectively. 
Primes denote the experimental reference pressure 
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of 2.000 atm. The compression k (Vv Vo)/V 
\lo \ ut pressure Pp: hi’ Is the corresponding com- 
pression nt 2? OOO atm. 1. @ (y’ lo) VY S is the 


area of the bore of the vessel at 1 atm: @ is the pres- 
sure coefficient of this area (1/8) (dS/dP); AL is the 
difference mn travel of the piston at pressure 7? in the 
measurements on the reference specimen and the 
unknown specimen calculated from 2,000 atm; Al’ is 
the difference of the positions of the piston at 2,000 
atm in the two measurements. The generalization is 
contained in the two summations, the first being 
earried out over all m components in the vessel 
except the pressure-transmitting liquid, and the 
second for all except the specimen and the liquid. 
All components of the compression specimen were 
weighed and the required volumes calculated from 
the weights. Known densities of steel and mercury 
were used, whereas that of the glass was determined 


by hydrostatic weighing. Volumes of the hydro- 


carbons were obtained from the weights bv use of 


specific volume-temperature data of Templin [7]. All 
these data were extrapolated to 21° (, and some 
probably not exceeding | ; eXists in 
The values used were, in 


uncertamts . 


the specific volumes used. 


cubic centimeters per gram, 1.070 for C,.H,.; 1.060 for 
The latter 


CoH; and 1.050 for all other specimens. 
Value agrees with a flotation density of n-C. H, 
obtained by Miller |S]. The error introduced through 
inaccuracies in these values ts believed to be less than 
the precision of the compression measurements 

\ll glass tubes used were prepared from a single 
length of tubing, and all were annealed simultane- 
The compression of glass used was that re- 
ported by (dams and Gibson | The compression 
of mereury was also taken from the data of Adams 
10} and represents only a slight modification of his 
The compression of steel did not 
volume of 


ously 


published data [6] 


enter the calculations because the steel 
was the same in all measurements. 

\bout | em, of each of the following hydrocarbons 
was available for study: n-CyHgs, n-CoHy. n- 
C,H, n-C. Hy, n-CxHx, and n-CyoH... The melting 


points, which were sharp, are given in table | 


ABLE |] Co 
( 
Call CoH ( 
yon' | { 
\ \ \ 
J is ‘ is . 
wis To wis s 
‘ 7 wil ; wun wil ‘ 
} | 1 | 
Ws 1 ite ry. ‘ 1 
‘ “ teu “4 * — { 14 
: a4 2 s } 
} ! lL / . 
a ’ ef aun (ane ware wan aan ss 
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3. Results and Discussions 


The results of the compression measurements on 
the hydrocarbons are given in table 1. Values of 
compression, AV Vy, referred to 2,000 atm are 
tabulated for each 1,000-atm interval; measurements 
at 1,000 atm appear as negative values because the 
data are based on the reference point at 2,000 atm. 
The data at 1,000 atm are less reliable than those at 
the higher pressures and are not considered in the 
following discussion. Two independent measure- 
ments, obtained at different times, were made on each 
specimen, and these data are tabulated as runs | 
and 2. The averages of the two measurements on 
each specimen are also tabulated. 

By inspection of the duplicate measurements 1 1s 
apparent that the errors are large. The increased 
error |11] is to be expected as the specimen has a 
volume of only about | em* out of a total of 25 em’ of 
material in the pressure vessel. Ordinarily a 10- to 
12-cm® specimen is used. In view of the decreased 
precision a statistical analysis of the data was car- 
ried out. This analysis corroborated a standard 
deviation of +0.0040 in—AV V5 and showed that, 
for a given pressure, no really significant difference 
exists between the compressions of the different 
hydrocarbons within experimental error. On this 
basis the data of table 1 (excluding those taken at 
1,000 atm) have been averaged, with the results given 
in the last column of the table. These average data 
are considered to be representative of all specimens 
studied To obtain compressibility data, the aver- 
age compressions of table 1 were fitted to the Tait 
equation [12] by means of least squares. The result- 
ne equation, 

AV /V,=0.1717 log, (2.500 2?) — 0.6272 2 
with /? in atmospheres, reproduces the average com- 
pression data within + 0.0005, which is well within 
the experimental error. The fit of eq 2 with the 
experimental data as averaged at each pressure 1s 


shown in figure 1. The corresponding values of the 
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KiGuRE 1] Plo ) la f equation for norma hudroca hon 


compressibility. at the various pressures are given in 


figure 2. 

It has been noted by Gibson |12] that the 
equation reproduces both volume and 
bility data on many materials with great accuracy 
wide range (see also {4]). Of im- 
mediate interest is the fact that this equation, 
although empirical, has been used very successfully 
in extrapolating from low to high pressures. The 
reverse process is used here, as data are desired at 


Tait 
compressi- 


over a pressure 


low pressures Krom the derivative of eq 2 it Is 
found that the compressibility is given by the 
equation 

dV O.1717 


VidP 2.303 (2.500-+-P) 


atm, the 
2G RW 10 


lf 2? < 2,500 volume compressibility is 
found to be atm Statistical analysis 
of the data indicates that the standard deviation in 
the value is 0.4 107° atm The true value may, of 
course, lie outside these limits 

lt is interesting to compare the present results 
with those of Miller |2], who measured the change in 
X-ray spacings wit lh pressure for the odd-carbon com- 
pounds n-C,H,y and n-CogH yo. In making this com- 
parison, it should be remembered that the even- and 
hvdrocarbons have different crystal 
structures at room temperature, so that close agree- 


odd-carbon 
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ment of the compressibility values is not necessarily 
to be expected. At room temperature the odd- 
carbon compounds from n-Cy,Hy to n-CygHyo are 
orthorhombic, whereas the even-carbon compounds 
from n-C\H,. through n-C.,H.o are triclinic; 
n-C gH. and n-CyoHy: are monoclinic [13]. Both the 
monoclinic and triclinic forms appear simultaneously 
in nm-C'o,H5, [1, 13]. Using Miiller’s finding [2] that 
the compressibility is essentially a two-dimensional 
effect? the linear compressibility values quoted in 
Miiller’s paper are equivalent to volume compressi- 
bilities of 21107° atm! for n-C.,H,y, and 5.510 
atm”! for n-CygHygo. These values represent the 
average compressibility at 500 atm; our value at the 
sume pressure for the even-carbon compounds is 
25.510°° atm”. 

The effects arising at chain ends lead to the expec- 
tation that the compressibility should 
slightly with increasing chain length in the normal 
hvdrocarbons. Bridgman’s data on the liquid hydro- 
carbons [3] up to n-CyH.. show this effect. The 
data of Cutler, Schiessler, and Webb [4] on n-C).Ho,, 
n-C) Hy, and n-C) Hy. in the liquid state also show 
a slight decrease of compressibility with length of 
chain. However, in the latter study it is clear that 
the effect decreases with increasing chain length and 
pressure, as would be expected. Bridgman |14] has 
reported data oa polvethylene of three different 
average molecular weights in the range 14,000 to 
38,000 that show no significant effect of molecular 
weight on compressibility, and verify the belief that 
the effect of chain ends on compressibility becomes 
negligible as the chain length increases im this par- 
tially amorphous and partially crystalline (mono- 
clinic) polymer. The fact that no effect of cham 
length on compressibility is shown by the present 


decrease 


lhe changes with pressure of the X-ray spacings found by M tller were stated 
to be almost entirely a result of the hydrocarbon chair moving closer together 
laterally 
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the crystalline hydrocarbons implies that 

1 an effect is smaller than the experimental error 
1 the Similarly, there does not ap- 
ear to be any large difference in the compressibili- 
of the triclinic and monoclinic forms of the 
ven-carbon hydrocarbons because the results for 
materials with 18 to 24 carbon atoms do not differ 
rreatly from for materials with 28 and 30 
carbon atoms per molecule. 


ti or) 
measurements 


those 


Mandel for advice 
statistical analyses of the 


The authors are indebted to .J 
in with the 
expel imental data. 


connection 


WASHINGTON, September 7, 
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Absolute Calibration of the National Bureau of Standards 
Photoneutron Standard: II. Absorption in Manganese 
Sulfate 


James De Juren' and Jack Chin 


An absolute method of calibrating a neutron source by determining the neutron ab- 
sorption rate in a manganese sulfate bath is described In this method slightly more than 
one liter of solution from the bath is irradiated in a strong neutron flux. The absolute 
activity of manganese in an aliquot of this activated solution is determined with a 42 flow 
proportional counter. Concurrently, one liter of this solution is added to the inactive bath 
and after stirring, the counting rate is measured with a dip counter As the number of 
radioactive manganese nuclei contained in the liter is obtained from the aliquot activity, the 
dip counter is calibrated for a known manganese activity. Then the bath is irradiated by 
the source and the dip counting rate at saturation vields the neutron capture rate by the 


manganese 


From a knowledge of the neutron capture cross sections and the manganese 


sulfate concentration, the fraction of neutrons captured by manganese is determined, As a 
result, the total neutron capture rate in the bath at saturation is obtained. A value of 
Q=1.25)* 10° neutrons per second was obtained for the National Bureau of Standards 
source (1954 The standard error is +2 percent 


1. Introduction 


As wide discrepancies between international neu- 
tron standards {1} * have been revealed by intercom- 
parisons, it was thought advisable to use two different 
methods for the calibration of the NBS standard 
neutron source. Both the emission rate of the NBS 
standard neutron source {2] and the thermal neutron 
density in a moderator [3] have already been deter- 
mined by methods in which the key measurement is 
the absolute counting rate of alpha particles from 
thin boron films. The purpose of the present work 
was to develop an independent method having an 
accuracy equal to or better than that obtained pre- 
viously. In the method selected the key measure- 
ment is the determination of the absolute activity of 
an aliquot from an activated manganese sulfate solu- 
tion, 

Two previous absolute calibrations of neutron 
sources have been made elsewhere in which manga- 
nese sulfate (MInSO), solutions were emploved to 
moderate and absorb the emitted neutrons. O'Neal 
and Sharff Goldhaber [4] used a simple method. 
The activ ity of a tank of Manganese sulfate solution 
is measured with a dip counter after irradiation to 
saturation with a neutron source and thorough stir- 
ring of the active solution. Powdered manganese is 
added to the solution and the mixture again trradi- 
ated to saturation After removal of the powder, 


the activity of the solution is measured and the 
absolute aAcUIVITN of the powder determined. The 
total activity of the powder is equal to (1—/)Q, 


where fis the ratio of the solution counting rates from 
irradiations with and without the powder, and @Q is 
the emission rate of the source When the irradiation 


Present lire West ise Electric Corp., APD, Pittsburgh, Pa 


Figures in bracket te the literature referem end of tl paper 
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is made with the powder present, f is the fraction of 
the emitted neutrons captured by the solution. 

The Swiss standard source has been calibrated in a 
similar manner [5]. Instead of using manganese 
powder, the activity of the solution is measured at 
two different concentrations to determine the ratio 
of the captures in manganese sulfate to water (hydro- 
gen). A solution of manganese sulfate is activated 
in an intense neutron flux from the D-D reaction 
The activity of the solution is measured, as before, 
with a dip counter. Part of the solution is  pre- 
cipitated as manganese carbonate and the absolute 
activity of a dried sample of this powder determined 
with a layer counter [6]. This measurement cali- 
brates the dip counter arrangement for this concen- 
tration in terms of the manganese activity present. 
Both methods require the determination of the ab- 
solute beta activity of a fairly thick laver of powder 
and the precision in the neutron emission 
claimed is of the order of 8 percent. 

In brief, the N BS standard Ra-Be \y,7) source con- 
sists of a beryllium sphere, 4 em in diameter, at the 
center of which a one curie capsule of RaBre is placed. 
As only gamma rays can enter the beryllium, the neu- 
tron emission rate is quite constant, unlike a Ra-Be 
an) source, Which has an initial growth of about 0.5 
percent a vear 


rates 


2. Experimental Method 
2.1. Outline of General Method 


In the present method the unknown source is sus- 
pended at the center of a cylindrical tank (approxi- 
mately 1 m in both diameter and depth) filled with 
manganese sulfate solution. After irradiation to 
saturation and thorough stirring, the activity of the 
solution is measured with a dip counter, Then in a 








ter experiment, slightly more than 1 liter of tank 
ution is irradiated in a strong neutron flux. <A 
mall aliquot (usually 0.075 ml) of this activated 
solution is pipetted onto a plastic film and its abso- 
lute activity determined with a 42 counter, One 
liter of this solution is added to the inactive solution 
the tank, the mixture stirred, and the activity 
measured with the dip counter. The number of 
radioactive manganese nuclei initially present in the 
tank is obtained from the absolute activity of the 
aliquot after correction for decay to the same arbi- 
trary reference time. Then the ratio of the satu- 
rated dip counting rate from the neutron source ir- 
radiation to the counting rate for the known manga- 
vields the saturated neutron counting 
This counting rate is converted 


nese acbuivIts 
rate Ih manganese 
to total neutron capture rate by a factor expressing 
the ratio of the total capture cross section per milli- 
meter im the bath to the manganese cross section per 
millimeter Let 

1) (, be the dip counting rate of the tank solution 
at saturation from the source irradiation, 

2) & 
activated solution is added to the tank, 

) Cy Be the absolute count ing rate of the 


quot of volume 0.075 ml 


ali- 


Then the ratio of the absolute manganese activity in 
the tank to the dip counting rate when the liter of 
activated solution is added to the tank and stirred is 


1000 , 
0.075’. C, 
I? may be interpreted as the absolute manganese ac- 
tivity in the tank per unit dip counting rate. The 


number of neutrons from the source captured per 


second by the manganese, (sin. IS 


If all the absorbing nuclei in the tank solution had 
capture cross sections varving inversely with the 
neutron velocity, the number of captures per second 


in the tank. 0. would be 


Nan Omn 7 %s Nuon 
(j= : i) val 
« \ «Mn: 

~¥axrn Ome 


where Vy, is the number of manganese (and sulfur 
nuclei per millimeter Vy is the number of hydrogen 
millimeter, and and om are the 
capture sections of manganese, sulfur, and 
hvdrogen for thermal neutrons. Neutron capture 
in oxygen is negligible for the photoneutron source 
\langanese is not a strict 1l/e absorber in the epl- 
small correction is 

Including the 


nucle: pet Tun, Fs, 


Cross 


thermal region and a ~1Q 
viven in the appendix 


term the equation is 


correction 


1.012 Nygo ona t+ Nagnos+ Nin 


ou 
. Ie 
é 1.012 Nuno: 


3b) 


s 


for the concentration used 


be the dip counting rate when the liter of 


A thin-walled Geiger counter, which had proved 
very reliable, was used as the dip counter. 


2.2. Description of the 4x7 Counter 


A 4x flow counter using a 90-percent argon, 10- 
percent carbon dioxide gas mixture was used for the 
absolute measurement of the specific activity. The 
counter has two hemicvlindrical collecting volumes 
separated by a plastic film 7.3 by 12.1 em. Stainless- 
steel wires, 0.0025 em in diameter and 2.5 em above 
and below the plastic film, extend parallel to the 
cylinder axes and are supported by kovar seals pro- 
jecting through the eyvlinder bases. One-quarter mil 
mylar films, coated on both sides with 5 to 10 ug em 
of aluminum, were mounted on a slide, which slipped 
into position between the hemicylinders (fig. 1). 

Pulses from the collecting wires were separately 
fed by cathode followers to Chase-Higinbotham non- 
overloading amplifiers operating at a gain of approxti- 
mately 30,000. The amplified pulses from each 
half of the chamber and the coincidence counting 
rate were recorded with a coincidence circuit having 
a resolving-time gate width of 2 ywsee. Several tests 
were made of the performance of the counter. An 
activated gold foil, 1 em? in area, was used to measure 
the relative counting efficiencies at different posi- 
tions on the mylar support. When the source was 
confined to a central region with boundaries 2.5 em 
from the bases and 0.8 em from the walls of the hemi- 
evlinders, the measured activities, corrected for 
decay, agreed within the standard deviation of 0.5 
percent. A calibrated Na®’ source (supplied by H. 
Seliger of the NBS Radioactivits Section) was 
measured with this 4% counter. The counting rate 
obtained, when corrected for decay and absorption 
in the mylar, was in agreement (within the experi- 
mental errors) with the original calibration. 

The counting rate as a function of collecting wire 
voltage and counting rate as a function of (integral) 
discriminator voltage characteristics of the chamber 
were Investigated with neutron irradiated manganese 
foils A collecting voltage of 1.950 v was selected in 
the center of a voltage plateau of more than 200 \ 
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The decay of the manganese foils was followed for 
several half-lives on six separate occasions. These 
measurements vielded a value of 223.7+0.2 min for 
the mean life, which is in good agreement with pub- 
lished values (7 | 

The 4x counting rate is equal to the sum of the 
individual counting rates from each half of the cham- 
ber minus the coincidence counting rate. Usually 
the collecting wires are connected together so that 
the 4x counting rate is obtained directly; however, 
knowledge of the individual counting rates gives 
more information regarding improper operation of 
the chamber. If a source is not thoroughly dried 
or the chamber not thoroughly flushed, the trouble 
is usually quickly detected by observation of the 
individual counting rates. 


2.3. Preparation of 4m Sources 


Considerable difficulty was initially experienced in 
preparing 42 sources, which, after drying, would 
adhere well to the mylar backing without flaking. 
The following technique proved the most satisfactory. 
An area 25 em? in extent was wet with a water solu- 
tion containing 1.5 percent of polyvinyl alcohol and 
0.2 percent of triton X-100. Triton X-—100 ts a 
potent spreading and wetting agent and the poly- 
vinvl alcohol an excellent binding material. An 
ultramicroburet, accurate to 1 part in a thousand, 
delivered 0.075 ml of the activated manganese sulfate 
solution onto the prepared area. The slide containing 
the manganese sulfate and a beaker of ammonium 
hvdroxide were placed under a bell jar for a few 
minutes until the manganese was precipitated. Then 
the beaker was removed and the bell jar evacuated 
through a cold trap to dry the film. The mass of one 
of these sources, determined by weighing, was 25.6 
mg, giving a value of about 1.0 mg/cm? for the 


average source thickness 


3. Experimental Data 


Five measurements were taken of the absolute 
counting rate of the 0.075-ml aliquots with a \4-mil 
mvlar backing, and one measurement was made with 
the aliquot sandwiched between two coated mylar 
films. For the backing alone a typical set of measure- 
ments, corrected for decay, was as follows: 

1) Single-chamber counting rate ~2,000 counts 
per minute 

2) Coincidence counting rate ~500 counts per 
minute. The top side of the chamber (facing the un- 
covered source) had a counting rate of about 6 per- 
cent greater than the bottom half. Table 1 gives the 
pertinent data for these six cases. 

\Ieasurements were made of the 4a counting rate 
asa function of mylar absorber thickness. To obtain 
av statistically significant effect, aluminum-coated 
sheets of mylar 7.5 mg/em* thick were added above 
and below the manganese sulfate source deposited on 
;-mil mylar film. The absorption data are shown in 
figure 2, A \-mil (~1 mg/cm’) mylar absorber causes 


TABLE 1. Calibration data 
——-—-- ~ —— 
A bso- 
in A bsorp- lute Dip 1.000 ¢ 
Source count tion count- Couns *R — = 
. ’ : ing rate, 0.0 < ¢ 
ing rate correc- ing rate, ( ‘ . 
tion Cur , 
Counts Counts! Counts 
On 4-mil mylar back min min min 
ing 2062. 8 1.014 | 2700.1 205. 7 1. 7X1 
Do 1.014 3068. 6 237.4 1. 724 
Do 3204. 4 1014 3340.5 252. 6 1. 763 
Do 3336. 4 1.014 3383. 6 257.2 1. 754 
Do 3170.9 1.014 3215.3 246. 6 1. 738 
In 44-mil mylar sand 
wich SO77. 6 1.020 3139.3 241.7 1. 732 
Average 1.74410 


* Each measurement of R, not including the error in the absorption factor 
has a statistical accuracy of about +0.9 percent, and the error in the average 
value of R is +0.4 percent rhe indicated uncertainties are expressed in 
terms of the standard error of the mean 


a reduction of (0.6+0.1%) in the counting rate. As 
the manganese sulfate source thickness is also about 
| mg/cm’, the correction for self-absorption will be of 
equal size; however, to allow for nonuniformity in the 
source thickness and variation in the area covered by 
different sources, the source self-absorption has been 
taken as (O.8+0.4%). All sources used should have 
self-absorption corrections falling within — these 
bounds. For the sources deposited on the mylar 
backing with no cover, the total absorption correction 
is (14+04%). At the concentration used and for 
the cross sections listed m the appendix, eq Sb 
becomes 
(J= 2.861 RC. (33¢) 
Using the average value of R from table 1, 
M—=4.990* 10°C. (4) 
With the manganese sulfate bath-dip counter system 


calibrated in this manner, the absolute calibration of 
the neutron emission rate from any source is obtained 
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from the dip counting rate at saturation. A small 
correction must be made for the self-absorption of 
thermal neutrons by the source when in the bath. 
When the NBS standard Ra-Be (y,n) source was 
placed in the bath, the saturated dip counting rate 


was 2.49, 1°) counts per second, giving a value 
(J=1.24, 10° neutron captures per second in the 
bath. In pure water the source self-absorption [2] 


umounted to 6,000+2,000 neutrons per second. 
For the manganese sulfate solution used, the slowing 
down length is about 4 percent longer than in water, 
and the mean life of a thermal neutron is only two- 
thirds of the life in water. As a result, the thermal 
neutron density in the vicinity of the source in the 
manganese sulfate solution is 60 percent of the value 
in water, reducing the source self-capture to 3,600 

Finally, the total emis- 
1.259 + 0.02;) « 10° 


1.500 neutrons per second. 
sion rate from the 
neutrons per second 


source is 


4. Discussion 


The major part of the final standard error in the 
calibration stems from the uncertainties in the cap- 
ture cross-section values for manganese and hydro- 
gen, Which cause a standard error of 1.2 percent in 
the capture factor. Recent measurements of the 
hvdrogen cross section by the pulsed neutron decay 
method yield a value of 0.333+0.003 barn [8], 
compared with the earlier value of 0.321+0.005 
barn {9}. The value of 0.333 barn, combined with 
the value of 0.329 +0.004 barn obtained by Hamer- 
mesh, Ringo, and Wexler [10] and 0.332 + 0.007 barn 
obtained by Harris et al. [11], would give a “best 
average” of 0.332 barn, compared with the 0.328 
barn used. Further work on the decay method is 
planned by von Dardel in order to resolve some of 
the uncertainties. However, an increase of 1.2 
percent in the hvdrogen cross section causes an 
increase of only 0.8 percent in the value of the neutron 
emission rate for the concentration of manganese 
sulfate used. 

A value of 13.2 barns was used for the thermal- 
neutron absorption cross section of manganese. As 
manganese has only one isotope, the activation cross 
section Is equal to the absorption cross section. The 
activation cross section of manganese was measured 
[12] by exposing a small cell of manganese sulfate 
solution in a known thermal neutron flux. After 
irradiation, the absolute activity of an aliquot of 
the solution was determined by te counting (with 
the same equipment used in the present experiment 
and Gece 13.1¢+0.39 barns was obtained In the 
second supplement of AECU-2040 (neutron cross- 
section compilation) an absorption cross section of 
13.2+0.4 listed for the thermal-neutron 
absorption cross section of manganese and 13.4+0.3 
barns for the activation cross section (the latter based 
ona gold absorption cross section of 98 barns The 
weighted average of the three values is accurate to 
0.1, barn 

Actually, systematic errors partially compensate 
other, when our own measured value of the 


barns Is 
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manganese activation cross sections is used. Any 
systematic error in the 4% counting technique, which 
causes an error in the value of Ci, causes a similar 
percentage error in the manganese cross-section 
result. These effects tend to oppose each other, and 
the net error in the source emission rate is only 
one-third of the systematic error. 

The emission rate of (1.25)+0.02;) < 10° neutrons 
per second obtained in this calibration agrees well 
with the value of (1.26;+0.03;) 10° neutrons per 
second obtained by a different method [2]. These 
two results give a weighted average of (1.25,;- 
0.02,;) 10° neutrons per second for the emission 
rate of the standard. 

The only American neutron source compared with 
a foreign source is the Los Alamos source. It was 
compared with the British standard at Harwell in 
1951 by three different methods [13]. As two of 
these comparisons differed by more than 4 percent, 
the exact ratio of the source strengths is not precisely 
known. In an effort to improve this situation, 
additional international comparisons are planned. 

A photoneutron source such as Sb-Be (y,n) could 
probably be calibrated absolutely to an accuracy 
within 1 percent by a modification of the present 
technique. As the gamma-ray energies are less than 
the photoneutron threshold of deuterium, a bath of 
manganese sulfate in heavy water could be used. 
Over 96 percent of the captures would be in manga- 
nese, and the error introduced by the cross-section 
values would be negligible. 


The authors express their appreciation to H. A 
Bright, NBS Analytical Chemistry Section, for the 
chemical analyses of manganous sulfate solutions; to 
H. H. Seliger for the use of his ultramicroburet and 
Na-22 source; and to R. Tilley and W. Bailey for 
assistance in performing the experiments. 


5. Appendix. Neutron Capture Above the 
Cadmium Cutoff Energy 


Manganese is not a strict 1/7 neutron absorber in 
the epithermal region; however, the resonances in 
the region above the cadmium cutoff energy are 
principally due to scattering [14], and the epithermal 
neutron absorption by manganese in solution during 
the slowing down is small. Walker has 
measured the neutron distribution from an Ra-Be 
av) source in water by means of manganese foils 
with and without 0.020-in. cadmium covers {15}. 


pre CeCSS 


For the space integrals (f° activity <r’dr) Walker's 
measurements give a ratio of (1.5 +0.1%) for the 
integral with cadmium covers to that with bare 


foils. A small BF, chamber gave a value of 0.71 
percent for the same ratio. As both boron and 
hvdrogen are good 1/r absorbers. 0.71 percent of the 
neutrons from the source are captured by the hvdro- 
gen in water the cadmium cutoff energy 
Manganese, therefore, has about twice the = epi- 
cadmium absorption of a 1/7 detector 


above 


The resonance escape integral, P,, may be used to 
evaluate the epicadmium captures. 


p “Bs Soa C, dE , 

, aia lal —ap (0) 
Eca >So; (’ tk 

where P, is the probability of a neutron escaping 


capture in slowing down from initial energy / to the 
cadmium cutoff energy, Lea; C;, ga; and os; are the 
concentration, absorption cross section, and seatter- 
ing cross section, respectively of the 7th element in 
the medium and € is the average logarithmic energy 
per neutron collision. In water and in man- 
ganese sulfate solution the slowing down from elastic 


loss 


collisions is almost entirely due to hydrogen, so the 
equation may be written 
Be Sal’, dE 
P,=exp- a (6) 
J Eca — H i’ 


with €=1 (for hydrogen). 
The integral is sensitive to the value used for the 
cadmium energy (usually taken as 0.3 to 0.5 ev). 


Bothe has derived an expression [16] for the proba- 
bility of absorption, a, of a neutron from an isotropic 
flux passing through a foil of thickness d and absorp- 
tion cross section o 


nod (nod 2k I nod . (7) 
where » is the number of nuclei per cubic centimeter. 
A graph of this function is given in figure 3 for 0.020- 
in. cadmium; where o,(/2)=o,(/) —7.0 barns is used 
for the absorption cross section of cadmium and 
a(f), the total section, is obtained from 
AECU-—2040. From the graph, £=0.5 ev appears 
to be the best choice for the cadmium cutoff energy. 
For water, letting 


Cross 


0 wu (iv) gon (th) (fe y,/F 
and 
"Eo gan (th) EY? ; 
Pp. exp | vali w - dk 
J Eca om =? 
2oan(th)F 
P,=exp- —_ : | ae 7 “a f (S) 
Osh ie! fer 
Using /y,=—0.025 ev, Ay~10° ev, oay(th) =0.328 
barn, and o.,;—20.0 barns, the value of P,—0.9927 


Thus 0.73 percent of the emitted neutrons are cap- 
tured by hydrogen in the epicadmium region. This 
figure is in agreement with Walker’s result for 
water. 

The concentration of the manganese sulfate used 
in the measurements was MnSQO,: 1.435 moles per 
liter, HO: 53.33 moles per liter 
The thermal neutron absorption cross sections used 
are orxn(th)=—13.2+0.2 barns, and o,s(th)=0.49 + 0.02 


cood 
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FIGURE 3 Probability of neutron absorption from an isotropic 


flux passing through YU.U20 in, of cadmium. 


barn. For this case, the escape probability is 
rae gk 7 Cee lk 
P. exp— Jan HT FaMn ‘Mn Tasl's ¢ (9a) 
J Eca onlu Ol 
*} O, dk 
Ps | exp = ad 
J Eca Fst 
"Eo gg, lh iw('s dE 
a Samn Mn ¢ i |le xp = Oa L ( : q (9b) 
JEca Ful | co Onn 


where &=1. 
Assuming that the thermal-capture cross sections 
vary as 1/#", the following values are obtained for 
the three factors comprising P,, respectively, 
P,=(0.9927) (0.99604) (0.99985). 

Manganese would capture 0.396 percent of the 
neutrons from the source in the epicadmium region, 
if it were a perfect 1/o absorber. From the man- 
ganese foil measurements, the manganese in solution 
will actually capture (1.5/0.72) «0.396% =0.825% of 
the emitted neutrons in the epicadmium region; 
hence 0.43 percent of the neutrons will undergo 
resonance captures in the manganese. From_ the 
thermal cross section data and the manganese sulfate 
concentre ation, 34.68 percent of the remaining neu- 
trons or 34.53 percent of the emitted neutrons are 
‘1/v”’ parce by manganese. Finally, (the sum of 
the 1/v and resonance captures) 34.96 percent of the 
emitted neutrons are captured by the manganese in 
the bath. For the concentration used, the manga- 
nese captures 1.2 percent more of the emitted neutrons 
than a perfect 1/r absorber, If this correction were 
neglected, the resulting error in the emission rate of 
the source would be 0.8 pereent 


315 








) 


Cy 


6. References 
J. ELughe Pile Neutron Research, p. 79 (Addisot 
Wesley Publishers Inc., Cambridge, Mass., 1953 
\. De Juren, D. W. Padgett, and L. F. Curtiss, J 
Research NBS 55, 63 (1955) RP2605 
\. De Juren and H. Rosenwasser, J. Research NBS 
52, 93 (1954) RP2477 
DD. O'Neal and G. Sharff Goldhaber, Phys. Rev. 69 
30S 1946 
Alder and P. Huber, Helv. Phys. Acta 22, 368 (1949 
\letzger, | Alder, and P. Huber, Helv. Phys. Acta 
21, 278 (1948 
Ht. Hopkins, Jr.. Phws. Rev. 77, 717 (1950); FE. Rab 
inowik Proc. Phys. Soc. (London) 63 [A], 1040 (1950u 
R. Schumann and A. Camilli, Phys. Rev. 84, 158 
1951 L.. M. Silver, Phys. Rev. 76, 589 (1949 
von Dardel and N. G. Sjostrand, Phys. Rev. 96. 1566 
1954 


316 


9] G. von Dardel and A. W. Waltner, Phys. Rev. 91, 1284 
1953 

10) B. Hamermesh, G. Ringo, and S. Wexler, Phys. Rev. 90, 
603 (1953 

11) S. P. Harris et al., Phys. Rev. 91, 125 (1953 

12) J. A. De Juren and J. Chin, Phys. Rev. 99, 1 (1955 

13) D. J. Littler, Ek. E. Lockett, A. D. MeRonald, P. W 
Mummery, and P. R. Tunnicliffe, British Atomi 
Energy Research Establishment Report, AE RE R/R 
776 (1951 

14] S. P. Harris, C. T. Hibdon, and C. O. Muelhause, Phys. 
Rev. 80, 1014 (1950 

IS! R L. Walker, I S. Atomic Energy Commission Docu- 
ment, MDDC 414 (Oct. 22, 1946). 

16] W. Bothe, Z. Physik 120, 437 (1943 


WASHINGTON, June 9, 


Vol , No. 6, December 195 


A New Bunsen-Type Calorimeter 


Ralph S. Jessup 


\ new Bunsen-tvpe calorimeter, using diphenyl ether as the calorimetric substance, ts 


deset bed 


Tests of this calorimeter show that it can be used to measure quantities of heat 


of the order of 38 calories with a precision of about 0.05 percent 


1. Introduction 


A new Bunsen-type [1]' calorimeter for use in 
measurements of small quantities of heat, such as 
heats of mixing of polymers and solvents is described. 
The calorimeter uses diphenyl ether (phenyl ether, 
phenoxvbenzene) as the calorimetric substance in- 
stead of water. The data reported in the literature 
indicate that Bunsen calorimeters using organic com- 
pounds as calorimetric substances have considerably 
higher sensitivity than that reported for the ice 
calorimeter [2,3]. The organic compounds that have 
been used for this purpose include acetic acid [4], 
anethole [5], diphenvimethane [6, 8, 9] naphthalene 
[10,11], phenol [7, 12], benzalacetone [13], and di- 
phenyl ether [14, 15, 16, 21]. The last-named sub- 
stance has the important advantages that it is quite 
stable, it can be easily prepared in a state of high 
purity, and its melting point (26.87° C) is conven- 
iently near room temperature 


2. Description of Calorimeter 


The calorimeter is shown schematically in figure | 
lt is similar in design to the ice calorimeter described 
by Ginnings, Douglas, and Ball [3], but has a larger 
central well to accommodate a “reaction vessel” for 
use in measurements of heats of reaction or heats 
of mixing. 

The central well of the calorimeter was made from 
2a plece of 1%, in. stainless-steel pipe that was turned 
down to a wall thickness of 0.02 in. at the upper end 
to reduce thermal leakage, and is hard soldered to a 
copper insert in the jacket cover, A, and to the brass 
piece, J, which constitutes the top of the calorimeter 
proper. The bottom of the well is hard- 
soldered to the stainless-steel pipe. The copper fins, 
FF, which are soldered to the outer wall of the well 
are partly to improve thermal contact between the 
well and the diphenyl ether, and partly for supporting 
the solid diphenyl ether, which is denser than the 
liquid. The copper spacers, (". between the fins at 
the lower end of the well are to promote uniformity 
in the vertieal distribution of heat transferred 
tween the well and the diphenyl ether outside of it. 
The outer shell, /, of the calorimeter is of stainless 
steel 4s in. in thickness and is soft-soldered to the 
brass piece, J The diphenyl ether occupies the 
space enclosed by the shell, /, the top, 7, and central 
well of the calorimeter, with the exception of the 


bra ss 


be- 


Figures in brackets indicate the literature references at the end of this payx 


space occupied by the pool, 7/7, of mereury in the 
bottom of £. 

An electric heater having a resistance of about 30 
ohms is wound on the outside of the cylindrical part 
of shell / of the calorimeter, and a heater of about 
10 ohms is wound on the hemispherical part of the 
bottom. Both heaters are cemented to the shell 
with glyptal lacquer, and are covered with aluminum 
foil to reduce radiation. They can be used sepa- 
rately or in series. The these heaters is 
described in section 4. 

When the calorimeter is in use the outer surface 
of the central well is covered by a mantle of solid 
diphenyl ether. Heat transferred between the cen- 
tral well and the diphenyl ether causes partial melting 
or freezing, with a consequent change in volume of the 
diphenyl ether and a flow of mercury into or out of 
the calorimeter through the platinum tube, W, and 
the stainless steel tube, J, so as to keep the total 
volume of liquid plus solid in / constant. The 
quantity of heat is measured by the volume or 
weight of mercury so transferred. If this is small, 
it can be determined from the displacement of the 
mercury meniscus in the 40-cm = graduated 
capillary, @, the steel valve, V7, being kept closed. 
If the volume of merceury is large it can be deter- 
mined from the combination of the change in weight 
of the contents of the beaker, PB, and the change 
in position of the meniscus in the capillary. 

The capillary, G, is horizontal so that the changes 
in the position of the mercury meniscus in it do not 
affect the pressure on the diphenyl ether in the 
calorimeter. Changes in barometric pressure might 
be large enough to affect the position of the meniscus 
significantly. To reduce any such effect the open 
end of capillary G is connected to a 1-liter flask 
immersed in oil in a Dewar flask. The pressure in 
this system is adjusted to atmospheric at the begin- 
ning of an experiment and has been shown by 
measurements with a closed-end mercury manometer 
to remain practically constant during the time of an 
experiment, 

The T-joint in the stainless-steel tube near the 
valve, V7, and the connections to the valve and the 
glass capillary, G, are made with Apiezon wax. 

The interchange of heat between the calorimeter 
and its environment is reduced to a minimum by 


use of 


glass 


>The calorimeter was originally constructed with the stainless-steel tube J, 
extending nearly to the bottom of the pool of mercury. Difficulties were en 
countered due to creeping of the diphenyl ether past the mercury seal between 
diphenyl] ether and the stainless-steel tube. This was overcome by replacing the 
lower end of the stainless-steel tube with the platinum tube, W’ [17] rhe nickel 
sleeve connecting tubes J and W is sealed to them with a low-melting-point glass 
made of equal parts by weight of soda glass, borax, and zinc oxide 
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surrounding the calorimeter with the aluminum 
jacket, 1, with its copper cover A. The jacket ts 
maintained at the temperature of the calorimeter 
within +0.001 deg C by means of thermoregulator T 
The temperature difference between the top of the 
calorimeter, 7, and the copper jacket cover, A, is 
measured by means of a 10-junction copper-con- 
stantin thermocouple (not shown in fig. 1) and a 
Diesselhorst potentiometer The coefficient of 
thermal leakage between calorimeter and jacket ts 
reduced by radiation shields (not shown in fig. 1 
of aluminum foil and by evacuating the 
between calorimeter and jacket to a 
about lu 

The 100-ohm electric heater, RP, the *>6-in. copper- 
nickel tube, ?, the evlindrical copper pieces, MZ, and 
My, and the conical copper piece, L, form a unit, 
the lower part of which fits close ‘ly inside the central 
well, ‘This unit, 
the calorimeter, 
similar unit having a reaction vessel instead of the 
heater, A, together with suitable tubes for introduc- 


space 


can be removed and replaced by a 


pressure of 


which is used only for calibrating 
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c drawing of calorime fer 


ing reactants. “Mhe purpose of copper pieces V4, and 
Mf, is to provide zones of good thermal contact with 
the wall of the well, and thereby trap heat that 
would otherwise escape by conduction along the 
leads, or by radiation and convection. The conical 
piece, L, fits into a conical hole in the jacket cover 
and forms a part of the jacket cover, and A/, forms 
part of the top of the calorimeter proper. 

The heater leads are brought out in such a manner 
that they are in intimate thermal contact with 
jacket ZL and calorimeter 14,. The potential leads 
are connected to the current leads midway between 
Land M As the total resistance of the current 
leads between Z and \/, is approximately 0.01 ohm, 
or O.OL percent of the heater resistance, the uncer- 
tainty in the heat developed between calorimeter and 
jacket is negligible 

The current through the heater was turned on and 
off by a double-pole double-throw switch similar to 
one deseribed by Osborne, Stimson, and Fiock [20] 
It is operated by a spring and a release that ts 
actuated by second signals from a standard clock 
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3. Filling the Calorimeter With Diphenyl 
Ether 


3.1. Diphenyl Ether 


The diphenyl ether used was purified under the 
direction of F. L. Howard by fractional freezing and 
distillation in the Engine Fuels Section of the Bureau 
The purity was determined by R. E. MeCoskey and 
G. T. Furukawa from measurement of the equilib- 
rium temperature of a solid-liquid 
various known values of the mass ratio of solid to 
liquid. The value reported for the purity is 99.9985 
mole percent and that for the triple-point tempera- 
ture is 300.03° K. It was later found [18] that the 
purity of this sample of diphenyl ether after heating 
to about 360° K in a tin-lined calorimeter had 
decreased to 99.9926 mole percent, 


3.2. Cleaning the Calorimeter 


Before shell fig. 1) was soldered to brass prece I, 
it was washed internally, and the central well and 
attached vanes were washed externally with organic 
After 
shell / was soldered in place, the interior surfaces 
were washed repeatedly with boiling distilled water 
until tests showed only a slight change in the re- 
sistivity of the water after such washings. The 
svstem was then heated internally and externally to 
un temperature somewhat lower than the melting 
point of the soft solder, and was evacuated for 
about a week at a pressure less than 107-* mm of 
mercury to remove water and = other adsorbed 
materials, 


solvents to remove traces of oil or grease. 


3.3. Introducing Diphenyl Ether and Mercury into 
the Calorimeter 


The method of introducing diphenyl ether and 
mercury into the calorimeter can best be explained 
by reference to figure 2. The liquid diphenyl ether 
was first put into a 1-liter spherical flask, FY and 
slowly frozen under vacuum. This process was 
repeated several times until most of the dissolved 
vas Was removed Then with the dipheny] ether 
frozen the flask was disconnected from the vacuum 
system and sealed to the 2-in. glass tube, 7° (fig. 2 
and connecting tubing, and connected to a vacuum 
pump and to the calorimeter as indicated. The 
system was then evacuated to a pressure less than 
10-* mm of mercury. 

During this procedure the end of glass capillary G 
fig. 1) was sealed, and beaker B (fig. 1) was replaced 
bv a short glass test tube containing mercury and 
closed at the top by a rubber stopper with a hole of 
the proper size to fit the glass capillary leading to 
valve Vo This tube was evacuated through a side 
tube above the mercury surface. Valve V was then 
opened so that the pressure in it and the glass capil- 
larv above it was reduced to that in the calorimeter 
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mixture for 








FIGURE 2 Schematic diagram to illustrate method of introduce 
ing diphenyl ethe nto calorimeter 

(less than 107' mm). Valve Vo was then closed, and 
air admitted to the glass tube, which replaced beaker 
B, thus filling the valve and connecting glass capil- 
lary completely with mercury. The beaker, B, con- 
taining a weighed amount of mercury, was then 
replaced in its normal position, as shown in figure 1 

An attempt was made to remove any residual gas 
in the diphenyl ether before introducing it into the 
calorimeter, as follows: 

With the vacuum pump operating, a small amount 
of diphenyl ether in flask / (fig. 2) was melted and 
allowed to run down into tube 7) where it was re- 
frozen by surrounding the tube with ice. The di- 
phenyl ether in 7’ was then remelted and allowed to 
run down into the U-tube at the lower end of 7. 
After this the U-tube was kept in an ice bath to pre- 
vent any premature flow of diphenyl ether into the 
calorimeter. The procedure of melting a small 
amount of the diphenyl ether in F, freezing it in 7, 
and finally remelting it and allowing it to flow to the 
bottom of 7, was followed repeatedly until the entire 
amount had been transferred from F into 7. 

Valve V (fig. 1) was then opened, and the desired 
amount of mercury was introduced into the calo- 
rimeter. After this, the temperature of the room 
was raised above the melting point of diphenyl 
ether, the entire amount of diphenyl ether in tube 
7 (fig. 2) and in the U-tube at the bottom was 
melted, and the liquid allowed to flow into the calo- 
rimeter. When the flow stopped, a small amount 
of liquid remained in the tube connecting to the 
calorimeter. Air was then admitted to the system 
through the line to the vacuum pump, raising the 
pressure to atmospheric. The connection to the 
calorimeter was then removed and the small opening 








into the calorimeter was closed by a conical stainless- 
steel pin forced into the opening by means of a screw- 
cap. Holes in ths sereweap permitted washing out 
anv diphenyl ether remaining outside the calorimeter 
proper in the space around the conical pin. 

It is believed that by the above procedure the 
diphenyl ether introduced into the calorimeter was 
substantially free from dissolved gases. Unfortu- 
nately, however, a small amount of air was later 
introduced into the system due to failure of the 
room thermostat during a weekend, resulting in an 
excessively low room temperature and consequent 
freezing of enough diphenyl ether to draw all the 
mercury from beaker B (fig. 1), together with a 
small amount of air, into the calorimeter. The part 
of the air that was trapped in valve V and tube J 
was easily removed, but that which actually got into 
the calorimeter could not be removed without a 
major operation, which it seemed desirable to avoid 
The air is normally dissolved in the liquid diphenyl] 
ether and apparently does not se ‘riously reduce the 
accuracy of measurements with the calorimeter, but 

t does cause considerable inconvenience in the use 
of the instrument 


4. Forming a Mantle of Solid Diphenyl 
Ether 


The inconvenience in the use of the calorimeter 
resulting from dissolved air in the diphenyl ether 
hes been greatest in attempts to form a satisfactory 
mantle of solid around the central wall. Because it 
is not possible to see into the calorimeter, some parts 
of the following discussion were inferred from the 
observed behavior of the calorimeter. 

Because diphenyl ether may undercool several 
degrees, it Is necessary to cool some part of the ma- 
terial in the calorimeter several degrees below its 
melting point to start freezing This is done most 
conveniently by pressing a piece of ice against the 
wall of the central well. If any considerable amount 
of diphenyl ether is frozen in this way, air in the gas 
phase accumulates in the calorimeter, as shown ‘by 
the fact that the change with pressure in the position 
of the mercury meniscus in the capillary, G, increases 
by a factor of as much as five in some cases. The 
excess sensitivity to pressure disappears due to solu- 
tion of air in _ remaining liquid, but mantles 
formed entirely by cooling with ice have never given 
satisfactorily consistent results in calibration experi- 
ments. It is inferred that this is due to air entrapped 
inside the solid in such a manner that it did not 
contribute appreciably to the apparent compressi- 
bility of the diphenyl ether, but that it dissolved 
the liquid when the adjacent solid is melted in a 
calibration experiment, thus causing the observed 
increase in volume of the system to be low 


In some causes the freezing of the diphenyl ether 


was started with ice, and the central well was then 
filled with water and freezing continued slowly by 
bubbling air through the water. This procedure 
did not always lead to satisfactory precision in the 
calorimetric experiments. In some cases it appeared 


that in this procedure the mantle did not spread 
uniformly over the central well, but tended to grow 
over only a small area, and finally to bridge over the 
space between the central well and shell / (fig. 1) of 
the calorimeter. 

The most satisfactory results were obtained with 
the following procedure: (1) A small amount of 
diphenyl ether was frozen by cooling with ice, and 
the calorimeter was left overnight to permit any 
gaseous air to redissolve. (2) The central well was 
filled with water, which was stirred by bubbling air 
through it; heat was supplied by means of theJdieaters 
wound on the outer shell / (fig. 1) of the calorimeter, 
and small pieces of ice were added periodically to 
the water in the central well. By proper adjust- 
ment of the heating current and the rate of addi- 
tion of ice, the rate of freezing of diphenyl ether 
could be made as low as desired, and the temperature 
of the water could be maintained from 0.5 to 2 deg 
below the freezing point. Under these conditions 
it Was expected that the mantle would spread over 
the central well so as to form a more or less uniform 
laver 

After the above procedure had been followed for 
about 1.5 hr, the addition of ice to the central 
well was stopped, and the current through the 
heaters was reduced gradually to zero as the tem- 
perature of the water in the well rose to the melting 
pot of the diphenyl ether Slow freezing of the 
diphe nvi ether was then continued over a period 
2 or 3 ds avs by bubbling air through the water in the 
central well, at first with the well filled with water 
up to the top of the calorimeter proper, and finally 
with about half of this amount in order to make the 
mantle thicker in the lower part where heat was to 
be added to it During this slow freezing a small 
current Was passed through the electric heaters on 
the outside of shall 

Even with a mantle formed in the manner tndi- 
cated, it was necessary to melt part of it by means 
of the heater, 722. before consistent results could be 
obtained. It seems probable that the part of the 
solid diphenyl ether which was frozen with ice con- 
tained entrapped gaseous air, and this part must be 
melted and the air dissolved in order to obtain 
satisfactory results 

It is believed that the difficulties encountered in 
forming a satisfactory mantle are due to the air in 
the svstem, and that if no air were present, a mantle 
could be formed much more quickly and easily 


5. Thermal Leakage 


The results of typical experiments to determine 
the thermal-leakage coefficient of the calorimeter are 
plotted in figure 3. The slope of the straight line 
corresponds to a thermal leakage of 0.0042 ) (0.0010 
cal) per minute for a temperature difference of 0.001 
deg C between calorimeter and jacket. That part 
of the thermal leakage due to conduction along me- 
tallie connections between calorimeter and jacket 
has been calculated to be approximately half of the 
observed total, so that the other half may be as- 
sumed to be due to radiation 
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TEMPERATURE DIFFERENCE, MILLIDEGREES 


ma eakage hetween calorime}lse and iachke 


FIGURE 3 The 


The observed thermal-leakage coefficient is a little 
larger than is desirable for measurements to a pre- 
cision of one or two hundredths of a calorie The 
“reaction period” for an electrical calibration exper- 
iment is about 30 min, and the temperature of the 
jacket varies in an irregular manner over a range of 
a few ten thousandths of a degree up to 0.001 deg 
It is found, however, that the observed rate of heat 
transfer remains more nearly constant over periods 
of 1 or 2 hr than would be expected on the basis of 
the curve in figure 3 and the observed variation in 
jacket temperature. This is undoubtedly due to 
damping of the effect of fluctuations in jacket tem- 
perature by the radiation shields between calori- 
meter and jacket, the relatively heavy calorimeter 
wall, and the lave of liquid between the calorimeter 
wall and the mantle of solid diphenyl ether. Al- 
though this damping of the effect of fluctuations in 
jacket temperature probably reduces the uncertainty 
in the thermal leakage, this uncertainty nevertheless 
is very probably one of the major sources of error in 
the calorimeter, particularly in experiments extend- 
ing over long periods of time 


6. Calibration of the Calorimeter 


The graduated capillary, @ (fig. 1), was calibrated 
by observing at various points along the capillary 
the lengths of 2- and S-em mercury threads with 
reference to the scale etched on the capillary , and by 
weighing the mercury removed from the capillary 
when the reading of the meniscus changed from 40 
em to zero. Corrections found for nonuniformity of 


321 


TABLE | Electrical calibration experiments 
) 1 atic 
Exper Calibration Deviation from mean 
me from mean 
ment factor Sow aoulon { LL exper 
l 4. 1435 0. 0021 0. 0012 
2 4.1452 O04 +. OOS 
} 4. 1458 +. OOO2 + OO11 
i 4.1473 +. O17 + 0026 
; $. 1433 Ou23 ool 
4. 1485 + (20 +, OOSS 
Mean 4+. 145 
tandard deviation of meat +0). (O08 
l 4. 1461 0. 0021 0.0014 
2 4. 1438 ‘LL TLD 
; +. 1418 W22 OU2Y 
) } i 4. 1421 oly 0026 
5 } : +. 1444 +. 0004 OO0S 
414M) + OO10 >. OOS 
| 7 4. 1443 +. OOO8 OOO4 
\ S +. 1449 + 0009 + OOO] 
Meat 4. 1440 
= ndard ad tion of ment oO 
4.1470 0. OOO 0. 0023 
y 4 +. 1445 + OO003 +. O00] 
; 1. 1436 OO009 OO 
i 41. 1445 + OOO08 + OOO 
4. 1425 OOO20 22 
I 4144 
st ial 1 I 0. OM 
Mean of all results 11447 
> 1 ia tor I LLL! 


the capillary bore did not exceed 0.2 mm (0.02 cal), 
and the mass of mercury required to fill the capillary 
was found to be 0.05239 . 0.000004 sdm) 
gem. This corresponds to an average capillary 
diameter of 0.71 mm. 

The calorimeter was calibrated electrically, using 
the 100-ohm heater, A, shown in figure 1, a Wolff- 
Diesselhorst potentiometer, a Wolff O0.l1-ohm stand- 
ard resistor, a 20000 20-ohm volt box, and a saturated 
cadmium standard cell maintained at a constant 
temperature of about 32.8° C by means of a temper- 
ature-control box [19]. The potentiometer, volt 
box, standard resistor, and = standard cell were 
calibrated in the Electricity and Electronics Division 
of the Bureau immediately following the calibration 
experiments reported, 

The results of three series of calibration experi- 
ments are given in table 1, expressed in absolute 
joules per centimeter displacement of the mercury 
meniscus in the graduated capillary, G. The results 
of each series were obtained on a separate mantle 
formed in the manner described in section 5. The 
amount of heat supplied to the calorimeter in each 
experiment was approximately 159 j (38 eal). 

Table 1 shows that the results of each series of 
experiments are in satisfactory agreement among 
themselves, the maximum deviation from the mean 
in any one series being about 0.07 percent, which 
corresponds to about 0.11 ) (0.02; eal) in the 159 j 
supplied to the calorimeter. However, the mean of 








the results of the first series differs from the means 
for the other two series by somewhat more than would 
be expected in view of the precision attained in the 
separate series. This suggests that perhaps the 
method of forming a mantle of solid does not entirely 
prevent the entrapment of air in the solid, and that 
he amount entrapped in different mantles may be 
significantly different, but that in any one mantle it 
s small enough and distributed uniformly enough 
through the solid to permit a relatively high precision 
it measurements with that mantle. It is to be 
expected, therefore, that to obtain a precision of the 
order of 0.1 j in the measurement of a small quantity 
of heat with the calorimeter it will be necessary to 
make a new calibration of the calorimeter for each 
new mantle 

If the individual results of all three series are 
averaged, the maximum deviation from the mean 
result so obtained of the mean for any one series is 
about 0.02 percent, corresponding to 0.03 5] 0.0leal), 
but the maximum deviation of individual results 
from this mean is now 0.09, percent (0.15 j, or 
0.04 cal). 

Combination of the calibration factor 4.1447 j/em 
with the value 0.05239 . gem for the mercury re- 
quired to fill the capillary vields the calibration 
factor 79.10 +0.010 (sdm) j gof mercury. This may 
be compared with the value 270.46 j/g¢ of mercury 
reported by Ginnings, Douglas, and Ball [3] for the 
calibration factor of the ice calorimeter. The 
sensitivity of the diphenyl ether calorimeter is thus 
about 3.4 times that of the ice caler'™meter. 

Combination of the value 79.10 j/g¢ of mercury for 
the calibration factor with the value 300.06° K for 
the freezing point of diphenyl ether under the 
pressure (114 em of mercury) in the calorimete: 
vields the value —A/// 7AV=35.6, bars per deg C 
This value is subject to some uncertainty because 
of the fact that the observed value of A/7/ AV may 
be affeeted by gaseous air in the solid diphenyl ether 

It is of interest to compare the calibration factor 
obtained in the present work with those reported by 


other investigators who have reported their results 


in terms of weight (or volume) of mercury expelled 
per calorie. The results on this basis are as follows: 


Author Calibration factor 
Hg/ca 
Sachse [14 0. 0364 
Holmberg [15 O488 
Klemm, Tilk, and Jacot It 052s 
(riguere, yn” ett ind Olmos [21 05200) 
Present work O5289 5 
7. References ~ 


[1] R. Bunsen, Pogg. Ann. 141, 1 (1870 

2) D.C. Ginnings and R. J. Corruecini, J. Research N BS 38, 
583 (1947) RP1796 

3] D. C. Ginnings, T. B. Douglas, and A. F. Ball, J. Re- 
search NBS 45, 23 (1950) RP2110 

1] L. kb. O. de Visser, Z. phvsik. Chem. 9, 767 (1892 

5] V. Grassi, Soc. Line. (5) 221, 494 (1913 

6] A. N. Shehukarev, I. P. Krivobalko, and L. A. Shehuka 
reva, Physik. Z. Sowjetunion 8, 722 (1934 

[7] A. N. Shehukarev, T. V. Ass, and N. L. Putilin, Méd 
exptl. (Ukraine), 1936, No. 6, 114 

8S} K.S Evstrop’ev, J. Phys. Chem USSR) 8, 130 (1936 

9) MI. M. Gordon, Zentralblatt 1938, IT 561 

10) C. C. Coffin, J. ¢ Devine, J. R. Dingle. J. H. Green- 
blatt, T. R. Ingraham, and S. Schrage, Can. J. Re- 
search |B] 28, 579 (1950 


11) A. Thomas, Faraday Soc. Trans. 47, 569 (1951 

12} S. J. Gregg, J. Chem. Soc. 1927, 1494 

13] Kk. J. Caule and C. C. Coffin, Can. J. Research [B] 28, 
639 (1950 

14] H. Sachse, Z. phvsik. Chem. | A] 143, 94 (1929 

15] Toivo Holmberg, Soc. Sei. Fennica, Commentationes 


Phvs.-Math. 9, No. 17 (1938 
16] W. Klemm, W. Tilk, and H. Jacobi, Z. anorg. allgem 


Chem. 207, 187 (1932 
17] R.S. Jessup, J. App. Phys. 23, 543 (1952 
IS} G. T. Furukawa, D. ( Ginnings, R. I \MeCoskey 


and R. A. Nelson, J. Research NBS 46, 195 (1951 
RP2191 

19] I FE. Mueller and H. F. Stimson, J. Researcl 
699 (1934) RP739 

20} N. S. Osborne, H. F. Stimson, and E. F. Fiock, Bs J 
Research 5, 429 (1930) RP209 

21) P. A. Giguére, B. G. Morissette, ar A. Wo Olme 
Car J. Chen 33. 657 L955 


NBs 13, 


WasHinecron, August 16, 1955 


322 








s 





Journal f Standards 


f Research of the National Bureau 


Vol. 55, No. 6, December 1955 Research Paper 2 


Mass Spectra of Thermal Degradation Products of Polymers 
Paul Bradt and Fred L. Mohler 


Polymers are degraded or evaporated from a tube furnace directly 
trometer, and mass spectra recorded as the temperature of the sample is increased 
ethylene degrades into paraffins, olefins, and diolefins and at 386 
Polyvinyl chloride degrades in two stages 

C 


tends to mass 684, Cy Hog 
it loses HCl and some benzene, and above 300 


into the mass spec- 
Poly- 
C the mass spectrum ex- 
From 127° to 300° ¢ 
it evolves a great variety of hydrocarbons 


The heavy fraction from rubber degradation evaporates over a wide temperature range, 
and molecules containing from 3 to 16 units of the monomer, C;Hs, are recorded rhe 


heaviest ion is CywHy; of mass 1,087 


ing from 1 to 9 units of the monomer, CyHg, the heaviest ion being C;.H 


Degradation of polyxylylene gives molecules contain- 


of mass 035 \ 


low polymer of polyphenyl was evaporated and molecules containing 2 to Ill monomer units 


were recorded 


1. Introduction 


The products of thermal degradation of polymers 
have been a subject of systematic research by 5. L. 
Madorsky [1, 2, 3]' and his colleagues in the NBS 
Polymer Structure Section. In Madorsky’s  re- 
searches the volatile products of thermal degrads- 
tion have been analyzed by mass spectrometer, and 
in most cases the average m¢ lecular weight of the 
nonvolatile degradation product has been evaluated 
by standard methods. In a few irstances, notably 
polymethyl methacrylate and polytetrafluoroethylene, 
the degradation product is almost pure monomer, 
and these techniques give a complete description of 
the degradation product. In most cases the non- 
volatile fraction is an appreciable part of the degra- 
dation product, and existing information on this is 
fragmentary and incomplete. 

In a previous publication by the authors [4] on 
polystyrene degradation, a technique is described 
for degrading and evaporating the products of degra- 
datioa directly into the ionization chamber of a mass 
spectrometer and recording mass spectra of all the 
products that are volatile at the temperature of degra- 
dation. The same technique has been applied to 
record the mass spectra of molecules of perfluoro- 
polyphenyl, which evaporates without degradation 
[5], 

The present paper gives a survey of mass spectra 
obtained when various other polymers are evaporated 
or degraded directly into the mass spectrometer. 


2. Experimental Method 


A small tube furnace extends to the entrance port 
in the ionization chamber of a mass spectrometer, 
and the sample to be studied is in a short length of 
capillary tubing supported in the tube furnace by a 
thermocouple, which measures the temperature of 
the sample. The amount of sample is of the order 
of 1 mg. The mass spectrometer is a Nier-type 60 
instrument of 6-in. radius of curvature, and the 
spectrum is recorded with a pen recorder by varying 
the magnetic field. The mass seale can be extended 


Figures in bracket ne ite the literature references at the end of this paper 


The heaviest ion was Cyl of mass 838 


as required by lowering the ion-accelerating voltage, 
which is normally 2,500 for molecular weights less 
than 300. The resolving power is about 1 in 300, 
which does not resolve unit mass in the high mass 
range encountered in these studies. Partial resolu- 
tion sufficient to identify individual ions of hydro- 
carbons extended to about mass 600. Ion currents 
are expressed in arbitrary units of scale divisions on 
the record at the maximum sensitivity used in each 
series of measurements. 

In each experiment, after obtaining a good vac- 
uum, the temperature was increased step by step, 
and at each step the mass spectrum was recorded. 
It required about 2 hr to cover the full range of 
spectrum at optimum recording speed. <A_ higher 
speed was often used for a quick survey to see 
whether there was anything to measure. Each ex- 
periment extended over several days, and after the 
sample was exhausted the ionization chamber was 
dismantled and thoroughly cleaned. 


3. Results 
3.1. Polymethylene 


A pure grade of polymethylene, made by decom- 
posing diazomethane in the presence of a catalyst, 
was supplied by the Polymer Structure Section. 
This polymer is similar to polyethylene, except that 
there should be no branched hydrocarbon chains. 
Madorsky and Straus [3] find that this degrades by 
breaking at random C—C bonds to give n-paraffins, 
l-n-olefins, and diolefins. About 97 percent of the 
degradation product is not volatile at room temper- 
ature and has an average molecular weight of 692, 
corresponding to about 50 carbon atoms. 

The sample of polvmethyvlene began to degrade 
into volatile compounds at about 335° C. At 352 
(' the hydrocarbon spectrum extended to mass 253, 
(\,H,;. The maximum peak was 179 scale divisions. 
At 369° C the ion currents were much larger but 
unsteady because of bubble formation. The mass 
spectrum extended to 460, Cy,Hg. As the sample 
became depleted the ion current became steady and 
at 386° C the heaviest ions appeared. Table 1 
gives some of the larger peaks in this spectrum, 
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LABLE |! Vass spectrum of de qg adation products of 
polumethulene at 386° € 
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which extended to mass 684, C,H... The maximum 
peak, at mass 57, was 103 scale divisions. 

The mass spectrum of table 1 shows a remarkably 
regular succession of peaks at intervals of 14 mass 
units from mass 253 to mass 684. For ions contain- 
ing less than 22 carbon atoms the larger peaks in 
each carbon group are, with a few exceptions, the 
alkvl ions C,H, In ions with more than 22 car- 
bon atoms the ions C,H become larger. The 
latter come presumably from the diolefin spectra 
The spectra gave no evidence as to the mass distri- 
bution in the degradation products except that there 
are molecules with 49 or 50 carbon atoms. 

The spectrum of table 1 was obtained from a 
sample that was nearly exhausted and this spectrum 
differs from the preceding spectrum at 369° © in 
the relative intensity of the pair of peaks C,H), 
and C,H,,4; in each carbon group. The 2n—2 peak 
is larger in the spectrum at 386° C, and the differ- 
ence becomes conspicuous al the heavs end of the 
spectrum This indicates that the polymethvlene 


becomes depleted in hvdrogen inh the process of 


degradation, 


3.2. Polyvinyl Chloride 


A sample of polyvinyl chloride, “Geon 101’, was 
supplied by the Bureau’s Plastics Section. Poly- 
vinvl chloride degrades Im a complicated manner 
At low temperatures evolution of HCL by the 


TABLE 2 Summary of mass spectra of de gradation prod icts 


ot pol yoirnul ( hloride 


lor urrents ile di 
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following reaction is predominant: 
H H H H 
( C—C—C »>—CH—CH—CH=—CH nC] 
= Cl BB Ci 


Some benzene is also evolved in this temperature 
range. At temperatures above 300° C the CH chain 
degrades to give a mixture of hydrocarbons. 

Table 2 lists the series of spectra obtained in one 
experiment in the range 127° to 389°C. It includes 
molecule ion peaks of HCl and benzene and the 
benzyl radical, and in the last column lists the 
heaviest ion recorded in’ each spectrum Peak 
heights are given in seale divisions 

The evolution of HCI began at about 127° C, and 
at 220° C the rate of evolution of gas reached the 
maximum amount for satisfactory measurement. 
The sample was held overnight near this tempera- 
ture, and then the temperature was slowly increased. 
Above 300 _* the evolution of HC} and benzene 
became small and a complicated hydrocarbon spec- 
trum appeared. The maximum peak in this spec- 
trum was the benzyl radical ion of mass 91. Above 
389° C the evolution of vapor decreased and a car- 
bonized residue was left in the sample tube. 

Table 3 gives the larger peaks in the mass spectra 
obtained at four temperatures. Below 300° C 
benzene is by far the most abundant hydrocarbon. 
The small peaks extending to mass 179 are heavier 
aromatics. The largest of these, at mass 128, is the 
molecule ion of naphthalene. The peak at mass 178 
may be the molecule ion of anthracene, but peaks at 
179 and 167 come from something else. At the 
highest temperature there are peaks at every mass 
number above 78, and bevond mass 200 all peaks 
become of the same order of magnitude. From 
mass 400 to mass 596 the spectrum of the partially 
resolved peaks appears as a very uniform sine wave 
with maxima at even mass numbers. The upper 
limit is set by resolving power, and the ion current 
is measurable bevond mass 596 
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TABLE 3 Partial mass spectrum of degradation products of 


poly in yl chlor ide 


lon currents in seale division 
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3.3. Crepe Rubber 


Crepe-rubber degradation has been studied by 
Straus and Madorsky [2]. The degradation products 
at 350° C are about 4 percent of monomer, 13 percent 
of dimer, and 883 percent of heavier compounds with 
un average molecular weight of 616. The monomer 
is isoprene, 2-methyvl-1,3-butadiene, of molecular 
weight 68, and the maximum peak is at mass 67. 
The dimer is mostly the evelie compound dipentene, 
and the maximum peak is at mass 68 [8]. 

Crepe rubber degraded in the range 236° to 353° C. 
The mass spectrum extended to mass 271 and by 
comparison with the known mass spectra of the 
monomer and dimer, the product appeared to be 
monomer, dimer, trimer, and tetramer in roughly 
comparable amounts. The experiment was not 
entirely because the ion currents were 
small. Failure to find the heavy molecules reported 
by Straus and Madorskyv [2] may 
Inadequate SeNSILIVIEN 

A sample of the heavy fraction of the degradation 
product with monomer and dimer removed was 
supplied by Straus. This was evaporated into the 
ass spectrometer and 25 mass spectra were recorded 
in the range SO° to 460° C. Table 4 gives some of 
the peaks in six of these mass spectra. These mass 
spectra are complicated, but, except for peaks at 
masses 69 and SI, the larger peaks come at intervals 


successful 


be ascribed to 


TARLE 4 Partial mass spectrum of the heavy fraction obtained 
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of 68 mass units, and have masses one unit less than 
the molecular weight of polymers of isoprene. 
Table 4 gives the abundances of these peaks relative 
to the mass 81 peak. The last row of table 4 gives 
the height of the 81 peak in seale divisions. ‘This is 
a measure of the rate of evaporation. 

At 80° C the vapor is mostly trimer. The peaks 
at masses 69, 81, and 135 are present but are rela- 
tively small in the dimer spectrum and must come 
primarily from ionization of the trimer. At 130° C 
the rate of evaporation is high, and molecules with 
$+, 5, 6, and 7 monomer units appear. After these 
evaporate, molecules with 8, 9, and 10 units appear. 
At 209° C there is rapid evaporation of molecules 
with 10 monomer units. These quickly evaporate 
and in the next spectrum record, molecules with 11, 
12, and 13 units appear. The heaviest ions were 
recorded at 264° CC. Peaks at 1,019 and 1,087 are 
barely visible and could possibly be accidental 
fluctuations. A comparison with the experiment on 
rubber degradation indicates that at 264° C there 
may be some thermal degradation of the sample, but 
at the lower temperatures this is negligible. The 
larger peaks in the lower molecular weight range 
must come almost entirely from ionization of the 
heavier molecules, 


3.4. Polyxylylene 


An experimental polymer, polyvparaxylylene, 
CH,-C,H,-CH,-CH,-C,H,-CH, . 

was supplied by the Polaroid Corporation. Experi- 
ments by Madorsky and Straus [6] showed that 
about 4 percent of the degradation product was vola- 
tile at room temperature and the heavy fraction had 
an average molecular weight of 661. Mass-spectro- 
meter analysis of the volatile fraction gave 78 mole 
percent of xylene, 7 percent of methylethylbenzene, 
1 percent of methylstyrene, 8 percent of toluene, and 
3 percent of benzene. 

When this polymer was degraded directly into the 
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PARLI ) Partial mass spectrum of de gradation products of 
pol yxrulyle ne at 310° C., 
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mass spectrometer it gave a complicated but repetl- 
tive Mass-spectrum pattern extending to mass 935. 
Table 5 gives the larger peaks in the spectrum ob- 
tained at 310° C. These peaks come at multiples of 
104 | mass units. That is, the larger peaks fall 
at odd mass numbers. The xvlene molecule ion of 
mass 106 is small, but xvlene amounts to only 3 per- 
cent of the degradation products 16]. The large 
peak at mass 105 must come from the ionization of 
heavier molecules. The complete spectrum includes 
other less abundant ions with 1 or 2 CH, radicals 
added to or subtracted from the chain of xvivlene 
units. Evidently the polymer breaks preferentially 
at the beta bond between the two CH, radicals, but 
it can also break at the alpha bonds. 


3.5. Polyphenyl 


The polyphenyvl was an experimental polyme 
synthesized by Max Hellmann of the NBS Polymer 
Structure Section from diiodobiphenyl heated = in 
contact with copper. The product was partially 
soluble in benzene, and the lighter soluble fraction 
was estimated to have a mean molecular weight of 
1.460, corresponding to 19 phenyl rings in each 
chain. The molecular weight of the insoluble frac- 
tion is unknown but is presumably much heavier 
than this 

When the insoluble fraction was slowly heated in 
the mass spectrometer to 500° C, only small impurity 
peaks were observed. The hydrogen mass range 
was not covered in any of these experiments. When 
the sample was removed it appeared to be completely 
carbonized. Apparently polyphenyl degrades by 
losing hydrogen atoms and not by breaking into 
smaller polypheny! chains 

The soluble fraction of polyvphenyl gave the results 
shown in table 6. Peak heights are given in scale 
divisions and many small peaks are omitted. The 
mass range below 230 was not covered in spectra 
obtained at 195° C and at 230° C. At the lowest 
temperature the sample gave off large amounts of 
unpolymerized biphenyl. (The ion of mass 163 is 
unidentified With increasing temperature heavier 
molecules appeared, but ion currents were always 


TABLE 6 Mass 8 pe ctra of the soluble fraction Oo; poly phen il 
lon current in seale divisions 
Mass lon 
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relatively small. At 195° C, molecules with 4 
phenyl rings are predominant, and at 230° ( 


molecules with 5 rings are most abundant, and the 
$ ring molecules are diminished. At 346° C and at 
390° C increasingly heavier molecules appear, but 
at the higher temperature rapid carbonization 
terminates the evaporation. 


4. Summary and Conclusions 


The experimental results with polvethylene are 
consistent with previous work on polvethvlene degra- 
dation, except that the upper mass limit 684 of 
table 1 is somewhat than expected. The 
average molecular weight of the heavy fraction had 
been estimated as 690 

Results on degradation of polyvinyl chloride 
should clarify existing uncertainty as to the mechan- 
ism of degradation |7]. There are two distinet 
stages in the degradation, with loss of HCI and 
benzene predominant at low temperatures and ran- 
dom breaking of the resulting CH chain occurring 
above 300° CC. The CH chain evidently degrades 
in a complicated manner above 300° C, taking up 
hvdregen atoms from the residue to form hydro- 
carbons with all saturation. There is 
no evidence of evaporation of chlorinated hydre- 


less 


degrees of 


carbons. 

The experiment on direct degradation of crepe 
rubber gave a range of molecular weights much 
smaller than expected, but this may be ascribed to 
The experiment on evapora- 
recovered from rubber 

of molecular weights 


inadequate sensitivity. 
tion of the heavy fraction 
degradation showed a range 
consistent with previous estimates. There are 
molecules ranging in molecular weight from 204 
3 monomer units) to 1,088 (16 units) and molecules 
with 9, 10, 11, and 12 monomer units are abundant. 
The average molecular weight [2] based on the 
melting point of a camphor solution was 616 or about 
9 units. 

Polyxvlylene seems to break at random = into 
fragments containing up to 9 monomer units. The 
range of molecular weights is entirely consistent with 
the average molecular weight, 661, of the nonvolatile 
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fraction. An interesting point, however, is that it 
has been assumed that polyxylylene actually degrades 
into monomer and dimer units and that these re- 
polymerize on condensation [6]. The present ex- 
periment gives no support to this theory, for degra- 
dation products are evaporated at pressures below 
10-* mm, and molecular interaction is negligible. 

The experiments with polyphenyl indicate that it 
does not degrade by breaking into shorter chains 
but loses hydrogen and carbonizes. A published 
research [5] on perfluoropolyphenyl shows that it 
too degrades by carbonizing. Low polymers with 
up to 11 rings evaporate without degradation, and 
the same upper limit is found for both the hydro- 
carbon and the fluorocarbon. 

The upper limits of molecular weights observed in 
mass spectra of four of the polymers are listed in 
table 7. These limits are presumably set by the 
consideration that heavier molecules will degrade 
at the evaporation temperature. A variable in- 
tensity factor also enters into consideration. The 
aromatic polymers are favored in this respect, as 
their molecule ions are relatively abundant. These 
are the heaviest hydrocarbons ever reported in mass 
spectra. The simple recurring pattern of these 
polymers leaves little doubt as to the molecular 
formulas of the heavy molecules. The polyviny! 
chloride product is an exception in that the degrada- 
tion product is hopelessly complicated. 

Three different types of experiments have been 
reported: (1) the direct evaporation of degradation 
products into the mass spectrometer, (2) the re- 
evaporation of degradation products recovered in a 
separate experiment (table 4), (3) the evaporation 
of low polymers without degradation (table 6). In 
the second and third types of experiments the com- 
plicated mixture is fractionated as the temperature 
is gradually increased, and this simplifies the re- 
sulting spectra. The results afford a basis for 


TABLE 7 [ pper mass limits of pol jmer mass spectra 


Polymer remperature Ion Mass 

°< 
Polyethylene 36) C eHes Os4 
Rubber 204 C oH: 1. ORT 
Polyxylylene 310 CoH: a3 
Polypheny! 300 Colts SBS 


observing the distribution of molecular weights in 
the material. Results are not quantitative, as there 
is no possibility of calibrating the sensitivity of the 
mass spectrometer in terms of pure chemical com- 
pounds. However, these experimental results are, 
in general, consistent with values of the average 
molecular weight of the degradation products as 
measured by the NBS Polymer Structure Section. 
There is no evidence that the condensed degradation 
products repolymerize to form larger molecules, 
but results on the degradation of rubber are in- 
conclusive on this point and deserve further study. 
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Heats of Combustion of Liquid n-Hexadecane, 1-Hexadec- 
ene, n-Decylbenzene, n-Decylcyclohexane, n-Decyl- 
cyclopentane, and the Variation of Heat of Combustion 


With Chain Length 


Frances Maron Fraser and Edward J. Prosen 


The heats of combustion of five highly purified long-chain liquid hydrocarbons were 
measured with a bomb calorimeter. The calorimetric data vield the following values for 
the heats of combustion, —A//c°(25° ©), of the liquid hydrocarbon with gaseous oxygen to 
form gaseous carbon dioxide and liquid water; | kilocalorie = 4.1840 kilojoule 


n-Deevlbenzene — 2340.58 + 0.42 keal mole 
1-Hexadecene — 2519.17 +0.44 
n-Hexadecane — 2557.15 + 0.42 
n-Deevievelohexane = 2497.90 +0.438 
n-Deevievelopentane = 2347.54 + 0.46 


It is shown that, in the 'iquid as well as in the gaseous state, the increment in the heat 
f combustion (and hence in the heat of formation) per CH, group added in the n-alky1 side 
chain is a constant except for the first two members of the series, and that this increment 
has the same value for each of these series of compounds 


1. Introduction measured. It was found that the measured heats 
agree, within the limits of experimental error, with 
This investigation is a continuation of the program — those calculated using the increment 156.26 per CH, 
a2 4 on the determination of some of the grroup added tn the n-alkvl side chain, 
thermodynamic properties of compounds of impor- This agreement shows that the heats of combustion 
tance in the national svnthetie rubber program. In of liquid hydrocarbons may be expressed by an 
the study of the processes of polymerization and de- | equation of the form A//,.=—A+Bm, m>1. <A is a 
polymerization, it is important to have thermody- | constant characteristic of the homologous series; B is 
namic data for the compounds In the liquid as well the constant for the increase in the heat of combus- 
as in the gaseous state, and data on the variation of | tion per CH», group added in the n-alkyl side chain, 
these properties with carbon chain length. Previous | and is the same for all these series; and m is the num- 
studies [5, 6) have shown that for several homologous | ber of carbon atoms in the side chain 
series of gaseous hydrocarbons the heat of combustion 


Increases by a constant amount upon the addition of ya Units of Energy and Molecular Weights 


each CH. group in the n-alkvl side chain: the first 
few members of each series are anomalous In an 
earlier paper [5] the variation, with chain length, in 
the heat of combustion was studied for the liquid 
normal paraffins The increment 156.26 keal per 
CH. group Was calculated. by the method of least 
squares from these experimental data This inere- 
ment, together with previously reported values for 
the heats of combustion of the l-alkenes /2, 22, 23, 24 
alkvibenzenes [19, 20) n-alkvlevelopentanes [21], and 
/ -alkvlevelol exanes |21) was used to calculate values 


for the heats ol combustion of higher members of 3. Method and Apparatus 


these series | ~ to the present time there has been 


The unit of energy upon which values reported in 
this paper are based is the absolute joule j , For 
conversion to the conventional thermochemical 
calorie, the relation, 1 cal—4.1840 }, Is used 

The molecular weight of carbon dioxide, the mass 
of which was used to determine the amount of reac- 
tion, was taken as 44.010, from the 1952 table of 
international atomic weights \7 


ho experimental check on the validity of this extra- The heats ol combustion of the compounds were 
polation Therefore, it seemed desirable to measure determined by burning the liquid samples In oxygen 
the heat of combustion of a higher member of each ina bomb, at constant volume. The procedure is 
series. In this investigation the heats of combustion | SUnilar to that used in previous investigations [9,10, 


ol l-hexadecene n-deevibenzene, n-deevlevelopen- 


| the atomic weight of carbon has recently been changed fror 
Lilie } -leey eve lohexane, and / -hexadecane were 9 12.011 [8 Che molecular weight of carbon dioxide was taker 
4.010 in th paper to retain cor tency with ilues previously reported | 
he Thermochemical Laboratory for the heats of formation of other hydrocarbor 
} | t ( | ofl svnthetic Rubber, u ind of carbon d cle Phe values of the heats of combustion per mole of hyd 
vith t (1 svnthnet hubbe Progra! rbon reported here can be ¢ rrected for the effect of the change to the new 
; , " the literature referet t the e1 pay tor eight of carbon by multiplication by the factor 44,011/44.01¢ 
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11,12). The samples were sealed in weighed glass 
ampoules to prevent evaporation before the com- 
bustion period. The filled ampoules were weighed 
and placed in a platinum crucible in the bomb, 1 ml 
of water was added to the bomb, and the bomb was 
flushed and filled with purified oxygen to 30-atm 
pressure at 25° C. The internal volume of the 
bomb was 385 ml. The sample was ignited by 
means of an iron wire fuse, 0.012 em in diameter and 
5 em long, coiled above the bulb. The amount of 
reaction was determined from the mass of carbon 
dioxide formed in the combustion. The heat of 
reaction was referred to the final calorimeter tem- 
perature of 28° C. Corrections were made for 
deviations from the standard calorimeter system, for 
the heat contributed by the ignition process of heat- 
ing and burning the iron wire fuse, and for the heat 
produced by the formation of nitric acid during the 
combustion. 

A series of calorimetric combustion experiments 
was made with each of the following substances: 
l-hexadecene, n-deevibenzene, n-hexadecane, n-dec- 
vievelohexane, n-decvlevclopentane, and benzoic 
acid (NBS Standard Sample 39g). All of the experi- 
ments were performed in as much the same manner 
as possible. The same calorimeter svstem, platinum 
resistance thermometer (NBS No. 485), and Mueller 
resistance bridge (NBS No. 24530) were used for all 
of the experiments. The experiments with standard 
benzoic acid were carried out to determine the 
energy equivalent of the calorimeter svstem. Some 
of these combustion experiments on benzoic acid 
were performed before and some after the experi- 
ments on the hydrocarbons, to allow for any errors 
which might result from changes in the calorimeter 
system with time. 

The heat of combustion of n-hexadecane had been 
previously determined [10] in the Thermochemical 
Laboratory. The measurement was repeated as a 
general check on the procedure. 

In several experiments, tests for completeness of 
combustion were made by analyzing samples of the 
product uses for carbon monoxide by the colori- 
metric method described by Shepherd [13]... In no 
case Was a significant amount of carbon monoxide 
found. A further check on the completeness of 
reaction was afforded by the ratio, 7, of carbon di- 
oxide produced in the combustion to the stoichio- 
metric amount of carbon dioxide calculated from the 
mass of sample used. The following average values 
of r were obtained in the experiments; the number in 
parenthesis is the number of experiments for which 
the value of r was determined, and the uncertainty 
is the standard deviation of the mean: 


O.0001] 
0.00008 
0.00005 
+ O.00012 
—~ 0.00009 
~ 0.00005 


0.99982 
0.99990 
1.00003 
1.OOO19 
1. QOO006 
1.00008 


l-hexadecene (5 
n-deevibenzene 5 
n-hexadecane (3 
n-deevlevelohexane } 
n-decvlevelopentane 5 
benzote acid (5 


4. Materials 


The hydrocarbons were from the API-NBS series 
of purified hydrocarbons. The assigned purities of 
these compounds in mole percent, as determined 
from measurements of freezing points, were: 1-hex- 
adecene, 99.84 +0.18: n-decvl benzene, 99 SS +0.10: 
n-hexadecane, 99.96 +0.04; n-deevlevclohexane, 
99.86 +0.11; n-deevlevclopentane 99.80 +0.18. 

The impurities in each compound are believed to 
be isomeric and present in such small amounts that 
their effect on the heat of combustion is less than the 
estimated uncertainty assigned. 


5. Experimental Results 


Table 1 lists the following data for the calibration 
experiments: The experiment number; the mass of 
the deviation from the standard 
the corrected increase in 
the cor- 


benzoic acid; Aeé 
calorimeter system; A/??,, 
temperature of the calorimeter system; q;, 
rection for ignition energy; Gn, the correction for the 
formation of nitric acid; and £,, the energy equiva- 
lent of the standard calorimeter system 


T ABLE | (a bration er pe ments 
Exe Muss 0 
mer be \ AR I 
No 
b h . hi 
51204 ] 0. JSODSO $1 D3u7u9 
2 1. 51221 Jit $4 ry LSYUS48. 4 
; 1. 5142 ' 2s47 $5.4 4 134S34 
{ 1. 5147 286751 sa | 2 13YS14.s 
51346 t Pa i | l 22 134814. 4 
t 1. 5315 ‘ INS $ ; 2 134820 
ri 1227 ; 2822 22.2 1ISYS15.8 
‘ 1. 5127 , JSS ; 22 SON]. ¢ 
Mi sUS21. 4 
Sta lard i } +4 
Exp I rm combustion experiment tt t 
bor 
Exper ry nbu 
rbons 


The standard calorimeter system was considered 
to be the actual calorimeter as it existed just prior to 
the start of the first calibration experiment except 
that the sample to be burned was not considered a 
part of the system. Thus Ae,, the deviation from 
the standard calorimeter, was computed from the 
heat capacity of the benzoic acid, 1.21) g° Cat 26.5° C 
minus the heat capacity, 0.03 )/° C ml, of that amount 
of oxygen at 30 atm and 28° C which was displaced by 
the benzoie acid) and the heat capacity of any parts 
of the calorimeter which were changed during the 
course of this investigation. Also included in Ae, 
are minor corrections for deviations of the average 
temperature of the experiment from the standard 
average temperature of 26.50° C (computed from 
the ealeulated change of heat capacity of the eal- 
orimeter with temperature) and for variations in the 
heat of the calibrating reaction due to deviations of 
the final temperature of the experiment from the 
standard final temperature of 28.00° C 
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Experi- 


Mass of carbon 


TABLE 2. 





Calorimetric combustion 


experiments 


Mean V slue of 
AUs with its 


’ (28° C q de a (2 
ment No dioxide Aes Af AS . . 4UB(2s standard devia 
tion 
n-Decylbenzent 
a ohm Ohms g CQ: gO 
l 9 S5874 11.8 0. 284165 39735. 7 31.5 41.9 13887. 0 
2 2. 84816 10.4 283114 3Y5S88. 3 35. 7 5.6 13885. | 
; 2 S925 10.6 285230) 39884, 3 $5.4 1.5 13886). 6 13886, 8S 
1 2 §7300 10.4 oN 4 30943. 5 35. 4 5.0 13889. 0 +0. 6 
2. 77388 10, 2 275735 38556. 5 t5. 2 2 13885, 3 
( > VI1OTE 10.¥ 280374 $0463. 8 35. 2 5.0 13887. 6 
1-Hexadecente 
l 2 67126 10.9 0). 2856083 30047, 7 35.8 6.0 14939. 0 
2 2 T2I7TRS 10.2 291792 J$OR01. 7 5.9 6.0 14941. 
7 2 66637 5 285222 S088 2. 8 36. 2 5.7 14942. 0 14941. 1 
1 2 65190 7.3 283676 30666. 0 35. 7 5.9 14941.4 0. ¢ 
5 2 (0973 7.3 27Y1 2 39029, 8 4. 4 5.7 14439. 4 
6 2. 6S7T06 7.3 287441 40192. 5 35.58 6.0 14942. 2 
n-Hexadecane 
l 2. 73431 5. 6 0. 206839 41506. 1 359 6. 2 15164.3 
2 2 67456 H8 20038 | 10603. 5 36, 1 6.3 15165. 5 
‘ 2 65107 6.1 INTTRO 10240. 4 35.2 5.7 15163. 5 15164. 5 
} 2 63905 5.1 286517 HOH. 7 36.1 6.2 15164. 7 +0. 4 
5 2 T7727 5.0 00919 42076. 4 44.58 6.4 15163. 8 
6 2 43637 5.0 244537 s6HUS8Y. 3 x6. 1 5.0 15165. 3 
n-Decylceyclohexane 
l 2.71781 4.1 0). 288249 10304. 6 35. 7 3.6 14815. 3 
2 2 32h 3.7 279230 30043. 4 35. 7 3 14815. 4 
2 59134 6 274887 38436. 1 5. 7 6.6 14816, 2 14814.7 
{ 2 60554 6.0 DSSS 5S {YUTU_S 35. ¥ 5. 6 14813. 1 +0. 5 
2 O57 RF 275319 N46, 5 35. 1 5.5 14813. 7 
n-Decyleyclopentane 
l 2. 48540 $4 0. 264257 $5049, 7 44.8 6.4 14850. 1 
2 2. 61842 3.5 TS461 5 oH. 2 7.1 14853. 4 
: 2 73536 1.9 *HG33 36, 2 23.5 14848. 3 14851. 2 
4 2. 73481 1.8 29148 35.8 23.8 14853. 9 +09 
2 AHB56 o.5 27271 $5.4 26. 2 14850), 4 
( 2 H4Hlt 6 28 1484 35. 5 | 24.6 14851. 1 


The energy equivalent, /,, was computed by: 


E.=((Q,m.44 Git Qn)/AR.|— Ae, 
where Q,—Al’g (28° C) is the heat of combustion 
of benzoic acid under the actual bomb conditions 
used in this investigation expressed in joules per 
gram of benzoic acid and m, is the mass of benzoic 
acid. The value used for the heat of combustion of 
standard benzoic acid under standard conditions at 
25° C is Qp (25° C) =26434.7 +2.2 j/g.° 
Under the actual bomb conditions and at 28° C. 
") =26431.8+2.2 j/g. 
I 


be | 


this value becomes Os, 28° ©) r 

Table 2 lists the following data for the calorimetric 
combustion experiments on the five compounds: The 
experiment number; the mass of carbon dioxide 
formed; Aéo, the deviation from the standard calori- 
meter system; AR,; Q (28° C)=(£,+ Ae) AR,, the 
isothermal heat liberated by the process which took 


J 


Chis value for benzoic acid, which is based on electrical calibration experi 
ments made in the Thermochemical Laboratory [11], was used to maintain con- 
sistency with values previously reported [10] for the lower hydrocarbons, which 
were based on electrical calibrations. This value differs slightly from the valuc 
certified for benzoic acid, Op (25° ¢ 26433.8+2.2 j/g [14, 15,. However, the 
uncertainty thus introduced is within the estimated over-all uncertainty. 





place in the bomb; qi; qn; and Al’g (28° C)=|Q 
(28° C)—qi—qn]/g COr, the heat of combustion in the 
bomb process expressed in joules per gram of carbon 
dioxide formed. 

In table 3 the following quantities are given for 
‘ach of the compounds: The mean value of —Al’, 
(28° C) expressed in joules per gram carbon dioxide 
from table 2, and in kilojoules per mole of hydro- 
carbon; —AL® (28° C), the decrement in internal 
energy in the ideal bomb process [10]; AHe° 
(28° C), the heat of combustion at constant pressure 
at 28° CC; —Alle® (25° C), the heat of combustion 
at constant pressure at 25° C; and A//,°, the heat 
of formation at constant pressure at 25° C. The 
uncertainties in the heats of combustion assigned 
to the values in table 3 were calculated using the 
relation, 


[(0.0083)?+ (20F,)? 
(0.01)7]}'/? 


uncertainty (2oAl,)” 


where o represents the standard deviation of the 


i mean for the specified measurement, and all de- 








TABLE 3 He als of comb istion and fo mation of the liquid h ydrocarbons at 25 + 
Ali 
‘ it \( it Zs” ¢ \( Zs~ ( AH f at 2 ( 
it 28°C it 2 ( 
,CO ki/ mole ’ é kj/ mole keal/ mole kj) mole keal/ mole 

ly Ltn 13886. 8 +2 O778. 5 1. 7 rt l l 791.3 l ‘ 4u743. 0 l t 340. AS 4), 42 219.12 +1.97 52.37 +0. 47 
1-Hexadecene 14041.1 2 6 10520. 9 1. 85 10518.0 +1.8 10538. 0 Ls 10540. 2 LS 2519. 17 0. 44 $29.41 +2. 00 78.73 +0. 

Hexadecane 15164 +9 1067S. 2 +1. 74 0675.5 +1. 74 166. 8 +1.74 10699. I 1.74 2557.15 +0. 42 $56. 35 200 109. OF +0. 48 

Decvieyelohexa 14814. 7 2 6 10431. 9 1.82 1429.0 +1.82 10449. 0 1.82 10451, 2 1. 82 2497. 90 +0. 43 418. 39 2 O68 100. 00 +0. 49 

Ly ke openta 14851 2.7 YSO4. 0 1.92 vSO1.3 +1. 92 20. 1 1. v2 YS22. 1 1.92 | 2347.54 +0. 46 368. 19 2. 12 88.00 +0. 51 


viations are expressed as percent of the measured 
quantity. The term 0.0083 percent is the assigned 
uncertainty in Q,; the term 0.01 percent is arbi- 
trarily assigned to allow for systematic errors in 
the measurements. 

The heats of combustion refer to the reaction: 
(b- 


C,H, (liq) ta)/4 O,(gas 


aCQ, (gas) +6/2 H,O (liq). 


The heats of formation refer to the reaction: 


aC (graphite) +6/2 H, (gas) =C,H, (liq 
The heats of formation were calculated from the 
heats of combustion using the values 68.317 


0.010 keal/mole [16] for the heat of formation of 
liquid water and 94.052 +0.011 keal/mole [14] 
for the heat of formation of carbon dioxide. 

The value obtained for the heat of combustion 
of n-hexadecane (liquid), AHe° (25° ©) =2557.1 
+ 0.42, is in good agreement with the value previ- 
ously reported [10] by the Thermochemical Labora- 
tory, AHe° (25° C)=2557.58 0.68 keal/mole. 
The agreement of this value with those of other 
investigators has been described in a previous publi- 
cation [5]. No previous measurements have been 


reported for 1-hexadecene, n-deevlbenzene,  n-de- 
evicyclohexane, or n-decylevclopentane. 
TABLI } ( om parison of ex perime ntal values of Alle 
Series 1-Alkeme Normal alk ylb 
Structure Vinyl-(CH H Phenyl-(CH 
Experimental Deviation | Exp : 
site liq calculated AN ! 
ul é 
) SI) US ] 
] 4.50 ] 
44. 40" +0). 1S 1) 10 1091.08 1 
SOM). OT 0. 26 Ol 1247.19 1 
} i. 72 34 12 1405. 4 0. 27 
1112. 87 0). 38 21 
Ww 234 5 0. 42 
14 2519. 17 44 7 
l 1&8 
R 16. 24 
At satu 


6. Heats of Combustion of Homologous Com- 
pounds 


The increment in the heat of combustion per CH, 
group added to the chain in the n-paraffins in the 
gaseous state has been shown [5] to have the constant 
value 157.44+0.05 keal/mole. This increment 1 
constant beginning with the third member of the 
propane. From a study of the heats of 
isomerization of branched-chain paraffin hydro- 
carbons, it has been shown that the average bond 
energy of a particular bond between two atoms is 
affected, within experimental error, only by the other 
bonds (or atoms) attached to these atoms. Bonds 
(or atoms) further removed have no measurable 
effect on this bond energy unless major steric effects 
are present (such as are revealed by the difficulty 
in making a Fisher-Hirschfelder model of the 
molecule). The constancy in the increment per 
CH, group for members of the series above propane 
is in accord with this observation, since the addition 
of a CH, group to propane to form n-butane forms 
the same kind of bond and neighboring bonds as do 
further CH, additions. In later reports [6,18,19] 
the increment of 157.44 was applied to other ho- 
mologous series of hvdrocarbons with n-alkvl side 
chains to obtain heats of combustion for highe1 
members of these series in the gaseous state. 

A constant increment per CH, group for numbers 
above propane was also found [5] for the n-paraffins 


Is 


Series, 


with calculated values using the equation Alle A Bm 


rem Normal alk yieyclopentan Normal alk yleyelohexanes 
i C yelopentyl-(C Hy) »—H Cyclohexyl-(C H2) »—H 
eviation wey tior leviallor 
n | Experimental . - nti experimental ane ntal 
rimental 8 ‘ veri i li, ‘ y ‘ ] 
calculated MI . wlculated MK - calculated 
1 nole ca nole vs mole alniw ¢ ws mote 
+2. 49 786.54 +0. 17 1.49 U6. SS OW 119 
0.23 941.14 +0. 1S O15 111 0.2 Oo. S0 
iW 1007. 50 +0). 22 os 1248, 23 ‘) iM) 
2 1253. 74 0. 2s i 1404. 44 0. 27 07 
Ol 1a). 7S ) ; ] 
$1 2347 { 4 ww 244 ") o4 1u 
78. 40 7S85. OF 35. 0Y 
24 Lt. 24 Lit. 24 
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The value of this increment is 
156.26 +0.05 keal/mole. Because of possible varia- 
tions in intermolecular forces, it was not evident 
that this increment could be applied to other ho- 
mologous series of hydrocarbons with n-alkyl side 
chains in the liquid state. The values obtained 
here for the heats of combustion of some higher 
members of such series, together with values previ- 
ously reported for the lower members in the liquid 
state, show that this value of 156.26 can be applied 
to these series as well. The constancy begins with 
the ethyl side chain. 

The constant B was re-evaluated, by the method 
of least squares, utilizing all of the available experi- 
mental data from the Bureau. The new value of 
B, 156.24, is in agreement with the value 156.26 
calculated from the n-paraffin series alone. 

In table 4 are given the experimental values for 
the heats of combustion of compounds of several 
series in the liquid state. These data were obtained 
from the following sources: the heats of combustion 
of benzene, toluene, ethylbenzene, n-propylbenzene, 
and n-butylbenzene in the liquid state [19, 20]; the 
heats of combustion of cyclopentane, methyleyclo- 
pentane, ethyleyclopentane, n-propyleylcopentane, 
cvclohexane, methyleyclohexane, ethyleyclohexane, 
n-propyleyclohexane, and n-butyleyclohexane in the 
liquid state [21]; the heats of combustion of 1- 
butene [2] and l-pentene [22] in the gaseous state 
together with heats of vaporization of 1-butene [23] 
and 1l-pentene [24]; the heats of combustion of 
l-hexene and I-heptene in the liquid state [25]; and 
the heats of combustion of 1-hexadecene, n-decyl- 
benzene, n-decyleyvclohexane, n-decyleyelopentane, 
and n-hexadecane in the liquid state from the present 
investigation. 

In table 4 are also given the deviations of these 
experimental data from the heats of combustion 
calculated using the equation Alle°=A+ Bm 
where A is the constant characteristic of the ho- 
mologous series, B is the constant 156.24 kcal/mole for 
each series, and m is the number of carbon atoms 
in the n-alkyl side chain. The values of A were 
calculated weighting the value of AH/c°, for each 
compound measured, inversely as the square of its 
assigned uncertainty. 


in the liquid state. 


| 
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Computation of Atomic Energy Levels: Spectrum of Singly- 
Ionized Tantalum (Ta nm) 


Richard E. Trees, W. F. Cahill, and P. Rabinowitz ' 


Energies and wave functions of atomic 


systems can be 
values and vectors of matrices in accordance with long-established procedures. 


characteristic 
( ‘oding for 


calculated as the 


carrying out these computations on the Standards Eastern Automatic Computer is described. 


Some preliminary values are 


tantalum (Ta 1). 


1. Introduction 


Calculations of atomic energy levels have been 
helpful in expediting the analysis of experimental 
data2* However, except in the simplest cases, the | 
numerical work is excessive for hand computation, 
and the calculations have had very limited applica- 
bility. It is possible to overcome this limitation by 
utilizing digital computers. A calculation of the 
levels of the low even configurations in the spectrum 
of singly-ionized tantalum (Tam) is given. This is 
a typical example of a calculation that is too difficult 
for hand computation, although fairly simple as 
compared with others in the field of complex spectra. 


2. Statement of Problem 


Spectra in which configuration interaction and 
spin-orbit coupling are both important are consid- 
ered; the basic theory for these spectra is well estab- 
lished.* ® According to this theory : energies of levels 
with the same total angular momentum (i. e., the 
same ‘J-value’’) are the eigenvalues of a matrix, 
which is here designated by (a;;). The J-value is 
omitted for brevity, but there will be one of these 
matrices for each J-value that is allowed, in the con- 
figurations of the spectrum under consideration. 
The configurations occurring in Taw are d‘, d 
and d? s?: these give rise to matrices of orders 9 > Q. 
1212, 2121, 1515, 15*& 15, 66, and 33 for 
the allowed J-values 0, 1, 2, 3. 4. 5, and 6, 
tively. 

Each matrix element is a linear combination of a 
limited number of ‘‘parameters’’, p®, multiplied by 


respec- 


known “coefficients” c’ (the parameters are the 
same for all J-values 
‘ L k f \ 
a .C3; P . (1) 
Present address: Weizmann Institute of Science, Rehovot, Israel. 
R. E. Trees and M. M. Harvey, J. Research NBS 49, 397 (1952) R P2378 
M. Gehatiah, Phys. Rev. $4, 618 (1954 


W. F. Meggers has recently located the c?2D(3/2) level of Hf at 30595, less 
han 100 K away from Gehatiah’s calculated position of 30692 

‘E. U. Condon and G. H Shortley, The theory of atomic spectra (Cambridgs 
University Pre ss, Londot 195] 

G. Racah, Puys os y 62, 4 8 (1942); 63, 367 (1943); 76, 1352 (1949). Referred 
to as R II, R III, ar aif vectively. 


given for the ene 





rgy levels of the spectrum of singly-ionized 


These linear forms can be derived in a straightfor- 
ward manner by use of methods developed by 
Racah (see footnote 5). For commonly occurring 
configurations, the calculation of any one element of 
the “coefficient matrix” (c) involves at most the 
multiplication of a few factors, which are given in 
readily available tables.° The parameter associated 
with a given coefficient matrix is defined in theory as 
a radial integral (Slater integral, spin-orbit integral, 
ete.) footnote 4). This can be evaluated if 
radial wave functions are determined by the method 
of self-consistent fields, but the work involved in this 
procedure is great. In the present calculations, the 
parameters are adjusted by least-squares to get best 
agreement with levels that have already been identi- 
fied experimentally. This limits the use of the theory 
to spectra in which the analysis of the experimental 
data is already well started. This limitation may 
also be removed when computers are used to carry 
out. self-consistent field calculations for complex 
spectra such as Ta 11. 

The coefficient matrices for the electrostatic inter- 
action in the configurations of Ta m are given in 
Rirand Rim. For each J-value there are 16 of these 
associated with parameters identified by the letters 
A, B,C, Gy, He, and @ (see table 2).7 In addition to 
these, there are three parameters ¢ associated with 
the coefficient matrices of the spin-orbit interaction. 
The latter coefficient matrices have been given for 
the configurations d*s* (see footnote 4) and for d' 
also.* The matrices for d* s have been computed by 
Trees; for the sake of brevity, the results are not 
given here. 

Recapitulating, there are 19 different parameters, 
and for each of the seven J-values there are 19 co- 
efficient matrices, one associated with each param- 
eter” T hough the orders of the coefficient matrices 


(see 


Bee: AuIsé th ere are often many elements, it will nevertheless require much time 
ind effort to set up these coefficient matrices correctly; probably automatic com- 
puters can be utilized in this phase of the work 

For a discussion of the parameter @ and other advances fn atomic theory not 
covered in footnotes 4 and 5, see G. Racah, Kungl. Fysiografiska S&llsk upets 
Handlingar N. F. 65, 31 (1955 Proceedings of the Rydberg Centennial Con 
ference on Atomic Spectroscopy). 

* H. Greyber, Ph. D. thesis, University of Pennsylvania, 1953 (unpublished). 

It is convenient to use the code of section 6 to calculate the g-values; to do this 
the LS-coupling g-values are included as an extra coefficient matrix, the ‘‘par- 
ameter’ associated with the g-value always being zero. Therefore, 20 parameter 
were actually used for Tal 
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are the same as those already given for the appro- 
priate J-value, there are very few nonzero elements 
in most of the coefficient matrices: the only coefficient 
matrices with large numbers of nonzero elements 
are those associated with the two parameters c(d' 
and ¢(d° s 


3. Eigenvalues and Eigenvectors of Real 
Symmetric Matrices 


Preliminary estimates of the parameters p“ are 
made by using methods given in section 7.1., and the 
matrix elements are then calculated according to 


formula (1). This can be done by hand computa- 
tion; or the coding described in section 5 can be 
used to generate the matrices on the National Bureau 
of Standards Eastern Automatic Computer (SEAC 
The eigenvalues and eigenvectors of the matrices are 
then to be determined. 

A SEAC code that will do this is based on the 
Duncan-Collar iteration procedure.’ The of 
this code with matrices up to order 1212 has been 
briefly reported."' The ean handle matrices 
with orders up to 2525. The major disadvantage 
of this code is that it requires too much machine 
time for regular use when all eigenvalues of several 
high-order matrices are desired Because of this, 
no effort was made to determine all the eigenvectors 
and eigenvalues for the matrices of Ta 1 set up with 
the preliminary values of the parameters. No 
matrices were set up with the final values obtained 
for the parameters. 


use 


code 


4. Least-Squares Adjustment of Parameters 


The matrices set up with preliminary values of the 
parameters p“ have the eigenvalues \°° and the 
eienvectors x determined the components of a 
vector will be indicated by deseribed in 
section 3. The eigenvalues are approximated as a 
function of the parameters by the following linear 
formula, which will be accurate to first order, in the 


is 


difference between the parameters and their pre- 
liminary values: 
( ae | i > 
x Quit Typ 2 
where 
| 
1m) — SY. pm, 
hemmed 
The matrix mod ifie t btain su 
WM. Kincaid, Quart. Apy Mat 5 
ID. I Mann, I Fano, W F. = I 
22, 764 (1954 
Appreciation is expressed to W. G. Hal ru SFA 
values in Ta u reported in section 7.2. Hea ht ; 
trum of tunesten (W 1) for G. Raeah: a ret t the latter work. w 
matrices with orders up to 95 & 25, has be . , On 
of his work, Hall has estimated that the 
eigenvalue and eigenvector about Smir ra . trix. & 12 
matrix, and 50 min for a 25 & 25 matrix 
kK. CGioldberg has a code that dets nes only 
ilar to that described by W. Giver ORNL 1574 j | . 
the eigenvalues of matrices up t le 0 Mt 
natrin e Al it 10 mina mied for a 2 ¥ 
being extended te ve the veetor t is expected that the exte 








The coefficients c; are the ones already used in (1), 
The linear formulas (2) can be set up by hand; or 
the coding described in section 6 can be used to 
generate them on SEAC. A least-squares adjust- 
ment can then be made to find new values of the 
parameters p“, which will make the linear forms (2) 
agree with a corresponding set of observed 
experimental values. 

A general-purpose orthonormalizing code has been 
set up for SEAC," and it can be used to carry out 
least-squares calculations. The coefficients dim, for 
fixed ‘“*k’’, are regarded as components of a vector in 
& space with the same dimensionality as the number 
of observations to be fitted. For successive k-values, 
each vector is orthogonalized to all preceding vectors 
with smaller k-values. The observed values, re- 
garded as an additional vector, are then expanded in 
terms of this set, and the residues that remain after 
the expansion will be the least-squares deviations. 
The process can be interrupted before the ortho- 
normalization is completed for all k-values, so that 
the least-squares calculation can be carried out, Say, 
first by omitting and then by including the last 
parameter, without making two separate calcula- 
tions. As many sets of observations may be fitted as 
desired in a single calculation. Consequently, differ- 
ent possibilities for correlating observations with 
theory may be tested. An observation may 
omitted from the calculation by inserting a new 
“weight vector” in which the component correspond- 
ing to the deleted observation is given the weight 
zero. The vectors must be “augmented” in these 
applications (see footnote 14), so that the least- 
squares values of the parameters will be recorded by 
SEAC (along with the deviations between theory 
and experiment This augmentation reduces the 
capacity of the code when used with the high-speed 
memory; with 18 adjustable parameters, a maximum 
of 22 observations could be fitted, whereas with 9 
adjustable parameters, 65 observations could be 
fitted. However, an alternative slightly slower code 
is available if external tapes are used, and in this the 
only limitation is that the sum of the number of 
observations and the number of parameters be less 
than 256. 

The Ta um problem was worked in an approxima- 
tion where first 8 and then 9 parameters were ad- 
justed (i. e., results were obtained with and without 
the final parameter @ included). There were 39 ob- 
but the number of these was reduced 
in successive calculations. As a check, both the 
observed values and the original deviations from 
these values were fitted. It required about 20 min 
of machine time for each of these calculations as com- 
pared to about 9 hr required in an earlier hand com- 
putation, where the normal equations were first 
obtained and then solved by elimination. The 
sults are described in section 7.3. 


best 


be 


sery ed values, 


P. Davis and P. Rabinowitz, J. Assoc. Comp. Mach, 1, 183 (1954), 
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5. Data Preparation—Matrix Inputs 


A code has been up that will generate all 
matrices with a single input of the parameter values. 
It the internal memory and can generate 
matrices of order up to 38X38. Checks in the code 
insure that the matrix generated is correct on wire 
and ready as an input for the coding described in 
section 3. For the Tau work, all matrices could 
be generated in 30 to 60 min of machine time. 

When a problem is first started, it is necessary to 
tvpe many coefficient matrices (20 for the Tau 
work) for each final matrix generated, and it would 
seem simpler to evaluate the elements (1) with a 
desk computer and type the final matrix directly. 
However, only the nonzero coefficients need 
tvped, and as there are usually very few of these, 
the total work involved is about the same.'® Once 
the coefficient matrices are typed, final matrices can 
be generated for as many sets of parameter values as 
desired. In any one problem, matrices must be set 
up with at least two sets of parameter values, so an 
appreciable amount of time will be saved." The 
sume coefficient matrices are used to generate the 
linear formulas (2) with the coding deseribed in 
section 6, and in relation to this work, the time 
required to type the coefficient matrices is a negligible 
factor. 

The coeflicient matrices used with this code must 
have no more than 56 nonzero elements, but larger 
coeflicient matrices can be inserted piecemeal by 
associating extra parameters (which will be given 
the same value) with each part. To allow for this 
possibilitv, the code has been written for 24 param- 
eters, Which is probably a larger number than would 
ever needed. The largest coefficient) matrix 
encountered ino the Tau problem = had* only 36 
nonzero elements. 


set 


uses 


be 


be 


6. Data Preparation Least-Squares Inputs 


A code has been set up that will generate the co- 
eflicients, formula (3), that enter into the linear 
formula for the eigenvalue, formula (2). Instrue- 
tions that determine the number of independent 
parameters desired in the linear formula outputs are 
first fed in; frequently, as here, parameters having 
the same symbol are regarded as identical in all 
configurations in which they occur, but other 
approximations are also used. All vectors for which 
linear formulas are desired are then fed into the 
internal memory; there is space for 582 vector com- 


required to set up the code 
, nd computers 


vcted that the same coding 


turally far in excess of the time 
Che expenditure of this time is 
will be used subse 


was tl 


juently in many similar prob! 


wilicient are used to 


Che last four digits of each ¢ specify the position in the 
rix; this leaves seven digits for the coetlicient itself. In the present work 
ly four significant figures were retained after the decimal point, because a 
result theoretically correct to within 1 Kayser was considered adequate with a 
mean deviation between theory and experiment more than a hundred times as 
‘ it 
rhe same set of coetti natrices can often be used to generate final matrices 
r other spectra that have the same, or related, configurations The coetlicient 
matrices used in the Ta 1 work also apply to the configurations d°, d? s, and d® s2, 
ind they have been used to generate the matrices for the first spectrum of ruthe- 
nium (Rud) as a start toward assisting with the analysis of this spectrum 


The 


for each 


ponents (all vectors of a 2424 matrix) 
coefficient matrices are then fed in, and 
one fed in (or for each group of consecutive matrices, 
depending on the instructions), as many coefficients 
of linear formulas print out as there are vectors in 
the vector input (the maximum number of possible 
outputs is 24). For the Tai work, the 39 linear 
formulas in nine parameters were computed in less 
than 60 min of machine time; originally this work 
| was carried out with a desk computer, and the time 
| required was of the order of 60 hr. 
| The coding just described groups coeflicients of 
| the same parameter only if they are in the same 
matrix. For use in the least-squares code of section 
t, coeflicients of the same parameter must be grouped 
together without regard to the J-value of the matrix 
in which they occur. The second part of this code 
out the preceding coeiflicient outputs, <A 
remaining uncorrected deficiency is that it does not 
give an output ready to be used as an input for the 
orthonormalization code of section 4; the augment- 
ing and check sums were inadvertently overlooked, 
so the results had to be outscribed, amplified, and 
reinseribed to get the correct input. This sorting 
can be carried out only for problems where the 
number of coefficients, formula (3), is less than 672. 
With 15 parameters (which may be the least number 
that can give good agreement in Tai) only 44 
observations can be handled (and 39 of SL possible 
levels have already been observed in Tai). This 
limitation is a very difficult one to overcome without 
altering the construction of the computer itself, 
However, it is always possible to combine two 
partial sortings and set up by hand the final input 
for the least-squares code. 

It would be convenient at this point to carry out a 
few manipulations with the sorted linear formulas, 
but the potentialities for doing this have not yet 
been fully explored. However, supplemental coding 
has been made that will evaluate the linear formulas, 
and also the differences between the formulas and a 
set of observed values, for as many sets of parameters 
desired, 


sorts 


as 


7. Results for Ta II 


7.1. Preliminary Estimates of Parameters 


The spectrum of Ta 1 is very complicated because 
of the presence of the three overlapping configura- 
tions of 5d', 5d° 6s, and 5d? 6s?. The analysis of the 
experimental data ” is far from complete, as only 
39 of S81 possible levels in these configurations have 
been identified. J-values are assigned to all identified 
levels, but the configuration and LS-coupling desig- 





§ With some amplification, the coding could be used to obtain nondiagonal 
elements of the transformed coefficient matrices (as well as the diagonal elements, 
formula (3)); it is estimated that matrices of order up to 18 X 18 might be trans- 
formed with this amplified code. Such a code would be useful in solving the 
problem referred to in footnote 6 because it is sometimes easiest to set up coeffi- 
cient matrices in a scheme that does not correspond to observed behavior, and 
then to use a simple transformation to arrive at a labeling that is more nearly 
ipproached in nature. 

*C. C. Kiess, G. R 


. Harrison, and W.J. Hitchcock, J 
(1950) R P2075. 
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ons have been made for only 26 of them. How- 
ver, the configuration assignments of six of the 
levels have been uncertain. Kiess (see foot- 


note 19) had assigned these levels to the °F and *P 
f the 5d configuration, but this assign- 

nt put the *P term below the *F term, whereas 

would be expected to be much higher on the basis 

a simple LS-coupling theory. It would have been 
easonable to assign the *F term to the 5d*® 6s con- 
curation, but this would leave the *F term in 
still unidentified and force one to conclude 
that the lowest levels of this spectrum, corresponding 
to the 5d? 6s? °F term, were completely unknown.” In 
any case, it was certain that configuration and LS- 
coupling assignuments would be approximations at 
best, so that it would be very hard to decide what 
experimental values should be correlated with 
theory in estimating preliminary values for the 
parameters, 

Because of this, the theory was first applied to the 
experimental values * in the simpler, almost com- 
pletely analyzed spectrum of singly-ionized hafnium 
(Hf 11), to get some idea of the values that would be 
expected for the parameters in spectra with 5d- 
electrons; no calculations had been published for 
these spectra at that time. The Hf 1 parameters 
were not evaluated by least-squares, but the values 
obtained agreed well with those recently published 
by Gehatiah (see footnote 3).” 

Preliminary parameter values were then estimated 
for Ta u by adjusting the Hf m values to get rough 
agreement with a few of the best identified levels of 
the Ta u spectrum, The magnitude of the read- 
justment was kept reasonably consistent with what 
would be expected from the known behavior of the 
parameters in spectra with 3d- and 4d-electrons. 
However, some of the parameters were poorly ad- 
justed, and acceptable preliminary values of the 


rims ¢ Os" 


jd Os 


In commenting on this, Kiess noted that too few strong transitions were | 
unidentified in his line list to give much weight to thi possibility privat 
communication rhe assignment made by Kiess was originally made also by 
Pt A. Klinkenberg, G. J. van den Berg, and J. C. van den Bosch, Physica 16, 
S61 (1950), but they later favored the second possibility that the lowest levels had 
not vet been identified (Physica 17, 167 (1951 theory was used to demonstrat 
the validity of the second alternative by den Berg Ph. D. The tn 
versity of Amsterdam, 1951 (unpublished His calculation is much 1 
fied to decide the point, however, and the calculation given here show 
ively that the first alternative is the correct one 

22 W. F. Meggers and B. F. Scribner, J. Research N BS 13, 625 (1934) RP’732 

22 Gehatiah’s parameter values for Hf u are 

B=435, C=1,530, G2=2,818, Hy =482, and ¢ =1,33¢ 

Analogy with results for 3d- and 4d-electrons suggests that the Ta u pat eters 
hould more nearly equal to the parameters of Ta tor W 1 (excepting the sy orbit 
parameter ¢ for the latter spectrum) than to those of Hf. Some recent unpul 
lished evaluations of parameters in these two first spectra are terest Racah 
has given the final parameter values that he obtained in the 5d‘ 6s? and 5d° 6s con- 
figurations of W 1 

B=371, C=1,900, G2= 2,700, Hy =434, ¢ = 2,089 and a= 46. 

His calculations are briefly described in footnote 7. The Tat calculations a 

for the 5d® 6s? configuration and neglect configuration interaction; they have 


been carried out with and without the Z (1+1) correction ee R. |} ree 
Phys. Rev. 92, 308 (1953) and also T. Kamei, Phys. Rev. 99, 789 (1955 
B=345, C =1,289, §=1,657, a=112 
B=278, C =2,100, [ =1,650 
rhe neglect of configuration interaction in the Tat calculations is probably not 


ount for the value of a being | 


nay ace irger than expected 





parameters were still obtained, so it is likely that 
this readjustment is not critical. Some details of 
the procedure will be given, however, mainly to 
bring out significant features of the parameter 
behavior. The preliminary parameter values are 
given in calculation I of table 1, and the final least- 
squares-adjusted values of the parameters are given 
in calculation IT. 

The preliminary values for parameters B and C 
were assumed to be 20 percent greater than the 
values in Hf 11, although an increase of only about 15 
percent would be expected from analogy with 3d- 
electron spectra. Calculation II shows not only 
that this overestimated B and C, but also what ts 
more important, that the final ratio of B to C in 
Ta um is about 8 percent smaller than the ratio in 
Hf i (it is 25 percent greater than the ratio for the 
W 1 parameters (see footnote 22). This ratio has a 
more constant value in spectra with 3d- and 4d- 
electrons (see footnote 22). The ratio may vary in 
spectra with 5d-electrons, but it is also possible that 
a nearly constant ratio would be obtained in these 
spectra if a higher-order approximation were used, 
as suggested at the end of this paper. 

The preliminary estimate G,—2,000 was made to 
favor the assumption that the lowest identified °F 
term belonged to the 5d*6s configuration, even 
though this assumption was not considered correct. 
Calculation II shows that G, should have been given 
the value found in Hf u. Analogy with 3d-electron 
spectra indicates that G@, should be slightly larger in 
Ta u than in Hf 1, but the increase is much less than 
expected for B and C. 

The value ¢=2,100 isa considerable overestimate; 
calculation IL indicates that a value slightly larger 
than the value applicable in Tat should have been 
used (see footnote 22), which is what would be ex- 
pected from analogy with 3d-electron spectra. 

To estimate the interaction parameter 7/2, the *F 
and *P matrices were set up in LS-coupling by using 
all the preceding preliminary parameter evaluations. 
The lowest terms of these matrices were expected to 
be very sensitive to the value of the parameter //s, 
so that errors in estimating the other parameters 
would not be too important. A value of 500 was 
assumed for //2, as it was expected that the parameter 
would have a value in Ta u that would nearly equal 
that in Hf ir. Analogy with the behavior shown in 
the 3d-electron spectra would lead one to expect that 
the value in Ta um might be smaller than in Hf 1, but 
as this integral depends strongly on the overlapping 
of d- and s-wave functions, the analogy would not be 
very reliable. The eigenvalues of the two 5x5 
matrices were obtained on SEAC w ith the coding de- 
scribed in section 3. The lowest *P term was found 
to be 800 K above the *F term, instead of 1800 Kk 
below as observed. It was decided that an extra 10- 
percent increase in the value of 7/7, would be justified, 
and that this, together with the inclusion of second- 
order spin-orbit effects, would explain the observa- 
tions with reasonable parameter values. 
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at raBLe 1. Configurations d? s*, d° s, and d* of Tait in intermediate coupling 
of 
Levels are identified by the J-values and by the label number, which is the order in energy sequence 
to Calculation I: Solutions of matrices for preliminary values of the parameters 
er Calculation II: Solution from linear formulas; parameters evaluated by least-squares, with the levels 18500 and 23295 omitted 
re . 
Calculation I Calculation II 
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a 25409 1. 155 
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[ HIST 1.33 6067 780) 1. 362 6224 37 
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K 7.2. Preliminary Calculation of Energy Levels given in calculation | of table 1. Partly because the 
' eigenvectors obtained in this calculation are based on 
()- : 


By using these preliminary values for the para- 
meters, matrices were set up for the Tar problem, as 
outlined in section 2. (The coding in section 5 was 
not then available.) The eigenvalues and eigen- 
vectors of these matrices were then calculated on 


SEAC. as described in section 3. The results are 


preliminary values for the parameters, no attempt 
has been made to give the conventional LS-coupling 
designation for the levels, and the levels have been 
identified by the J-values alone; the assignments 
given to some of these levels bv Kiess (see footnote 
19) are, however, roughly confirmed. 
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the preliminary parameters were esti- 
any effort to assign the 
levels to apy theoretical 
ood correlati m can be established between 
levels and calculated eigenvalues, 


corresponding observed and cal- 


ul making 
thre observed 


herey 


veen the 


4] val ICS There ure, however, three observed 

with J=2 (16424, 17168, and 18500) in a 

ion where there are only two calculated (16576. 
7966): the g-values indicate that the level which can- 
be correlated with theory is the one observed at 
SO00 Similarly, there are three observed levels 
with J 25295, 23620, and 24870) in a region 


where two are calculated (23976 and 24563): in this 
case the g-values indicate that the level which cannot 
be correlated with theorv is the one observed at 
23295. The levels with J—3 observed at 12435 and 
14581 were correlated with the respective calculated 
values 13679 and 13322. This is the only instance 
where the calculated and observed values were not 
correlated in the same with respect to 
energy. This correlation made when it was 
found with better parameter estimates that the linear 
formulas for the eigenvalues indicated that thei 
energies would cross. This correlation could have 
been made more easily by noting that it makes the 
g-values agree much better, but this criterion would 
also call for reversing the correlations of the fifth and 
sixth levels with J=2, and this was not done. With- 
in the limits of accuracy of this preliminary calcula- 
tion, the only other correlations that seem uncertain 
are the ones connecting the level observed at 12705 
with the level caleulated at 12375, rather than the 
ones calculated at 11206 or 14675, and also the cor- 
relation of 17231 with 17597 rather than with 18354 
Both these correlations worked well in subsequent 
least-squares calculations However, the 
squares calculations show that agreement would be 
unproved if the level observed at 16424 were disre- 
garded, rather than the one at 18500. It is likely that 
other possibilities still remain, that might give better 


sequence 
was 


least- 


agreement 
7.3. Least-Squares Estimates of Parameters 


Linear formulas for the eigenvalues were set up by 
hand (the coding of section 6 was not then available 
in an approximation with nine adjustable parameters ; 
the parameter a, which was neglected for the pre- 
liminary estimates, is the additional parameter 
Least-squares calculations were then carried out on 
SEAC, with and without the parameter a, as de- 
scribed in section 4. These calculations will be 
described briefly; the parameter values that were 
obtained are yiven in successive columns of table 2, 
labeled to correspond to the descriptions. 

(1) The parameters were determined by using all 
39 observed levels with the correspondence between 
theory and observation the same as that already 
given in calculation | of table 1. The purpose was to 
verify, what was already fairly certain, that no 
readjustment of the parameters would give a good fit 


for all 39 observations In this calculation the level 


Least-squares evaluations of the 


@ 32, d3 s, and d* of Tar 


TABLE 2 parameters of 


Che parameters A determine the centers of gravity of the three configurations 


C, ¢, and @ would be specified independently in all three configurations in 
higher approximations than the one given her Che interaction parameter 
f d‘—d? s?), has heen taken equal to the exchange integral Ge(d@s ind desig 
nated G2; the two interaction integrals I2(d*—d s) and Iie(d° s—d? s2) are taken 


equal and designated // 
rhe first of a pair of determinations grouped together omits the L(1+1) cor 


rection, the second includes it A correspondingly numbered section of the text 
specifies the observations included in making each pair of calculation 
Calculation 
I ram t 
| I] Ill I\ 
1(d'? s 1am) 18303 20272 20168 2QOTSS 20140 A446 YUSSS 
lid? s 12843 12892) 13744 13785) «13614 13641 13636 13057 
lid MS MOS M777 HAUS HOS H41 2HS30) 0% 
/ 11 124 is} 2 74 tert i7v Hy 
( 191 2USS Is41 1945 187 1947 S24 212 
(; 24H 2476 2u72 Sis Jeod Pee ELT) 2024 2UIS 
é 1US0) 20S 1776 S00) 1786 17S Is43 Is74 
Il 107 MI2 38 ‘41 538 5Y v2] ‘44 
a 35 1u Ww 
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at 18500 had a deviation of 3817. and the deviation 
for the level at 23295 was 1294. The level at 
12831 had the next largest deviation of 1316, which 
is an order of magnitude smaller. The mean devi- 
ation was +1048; the mean deviation, exclusive of 
the two badly fitting levels, was +515. 

The inclusion of the L(L+1 
a negligible improvement in the mean deviations, 
and the correction was negative (a 35), wherees 
all the evidence so far accumulated indicates that it 
should be positive. With the L(1- 1) correction the 
mean deviation for all levels was reduced to 1.038, 
and the deviation, excluding the two badly fitting 
levels, was 510. 

Il) The levels at 18500 and 
from the least-squares calculation: the deviations for 
this calculation are given in calculation Il of table 1 
The mean deviation of this calculation is 371: As 
this is considerably than the value 515 for 
calculation (1), deseribed in the previous paragraph, 
it indicates that the parameters of the latter caleu- 
lation were influenced by the inclusion of the two 
badly fitting levels. Comparison of the parameters 
themselves also indicates this In this caleulation 
there are three deviations for the levels 16424, 
I4581, and 12831) that exceed twice the mean devi- 
ation of all levels fitted. 

Again, when the L(L+1 
the mean deviation is reduced only slightly (to the 
value and the correction Is negative (a 19) 
Most of the improvement produced by including the 
L(L+1) correction can be attributed to the slightly 
better agreement obtained for the badly fitting level 
12831, which has its deviation reduced to 858; if the 
mean deviation is computed without this level the 
value +371 would be replaced by ~ 341, and the 
value +364 replaced by 340. The inelusion of 
the L(L+1) correction is therefore not justified. 

As calculations II] and IV indicate, the inclusion 
of three badly fitting levels has not had much effect 
on the agreement obtained for the other levels. 
Also. there seems to be no possibility of including 


correct ion produced 


23295 were omitted 


less 


correction is included, 


364) 


340 


of 


Ss in 
ter, 


COT 


ext 


yn 
at 


the Z(L+1) correction and getting an appreciable 
improvement. The best calculation that can _ be 
offered at this time is that in calculation II of table 1. 
Possibilities for improving this calculation will be 
given in the discussion. 

(111) The levels 18500, 23295, and 12831 were 
omitted from the least-squares calculation. The 
mean deviation was +335, which is only a little less 
than the value +341 obtained in (II Inclusion 
of the L(L+1) correction produced negligible im- 
provement, as the deviation was reduced only to 
- 333; the correction was still negative (a 10). 
It would have been difficult to justify the omission 
of the level 12831 if this caleulation had given much 
better agreement. 

(LV) Eleven levels were omitted in this calculation. 
These omissions consist of the two pairs, 18500 and 
22928, and 12435 and 14581, omitted partly because 
the calculated eigenvalues were close and it seemed 
possible, therefore, that the linear formulas might 
not apply; the two triplet groups 23295, 23620, and 
24870, and 16424, 17168, and 18500, which were 
omitted because of uncertainties in the correspond- 
ence between theory and experiment; and the levels 
23381 and 12831, omitted because they had rather 
large deviations that tended to be improved by a 
negative L(L+1) correction. The mean deviation 
of this calculation was 
the same terms in calculation (II) is +200). The 
only large deviation was —491 for the level 14628. 
With the ZL(L+1) correction included, the mean 
deviation was reduced to 151, due mainly to a 
decrease in the deviation for the level 14628 to 

306. Again the sign of the correction was nega- 


tive (a 35 
8. Discussion 


Results obtained in spectra with 3d- and 4d-elee- 
trons indicate that when the L(L+1) correction is 
used the mean deviation between theory and experi- 
ment will generally be less than 200 K, and that the 
sign of the correction will be positive.“ Calculation 
(11) of table 1 should therefore be improved because 


Racah (see footnote reports a mean error of 411 without the L(1+1) cor 

m, and an error 229 with the L(L+1) correction in W 1; as indicated in 

itnote 22, the sign of the L(1+1) correction is positive By coincidence, the 

ean error of calculation II, table 1 is 412 K, nearly identical with the corre 

| ling W alt howe no appreciable improvement was obtained by the 

use of the L(L+1).correctior In Hf u, Gehatiah (see footnote 3) obtains a mean 
rror of 257 without the u f this correction 


170 (the mean deviation of 


a better significant agreement between theory and 
experiment can probably be obtained. 

It is fortunate that the experimental analysis of 
the Ta 1 spectrum will be extended by members of 
the NBS Spectroscopy Section. The calculation 
given is already accurate enough to furnish con- 
siderable help. This analysis may, in turn, help 
solve the problem of the two excess levels pointed 
out previously. Ina preliminary survey, the switch- 
ing of J-values has been ruled out, as none of the 
levels concerned are poorly identified, The identi- 
fication of new levels will help to eliminate some of 
the arbitrariness involved in correlating the observed 
and calculated values. With more observations a 
calculation with extra adjustable parameters (sug- 
gested below) will also be better justified 

Eigenvalues and eigenvectors of the matrices set 
up with least-squares values of the parameters will 
be obtained on SEAC after the coding now under- 
way is completed (see footnote 13). The g-values 
obtained should be much more accurate than those 
already calculated, and if this is so, the uncertainties 
in making correlations may be resolved. Knowl- 
edge of the exact eigenvalues will show whether or 
not any of the large errors of the calculation are due 
to limitations on the validity of linear formulas for 
the eigenvalues. The linear formulas obtained from 
these new solutions should be accurate enough so 
that eigenvalues can be redetermined from them for 
any additional small changes that are made in the 
parameters, without setting up new matrices. 

The agreement between theory and experiment 
may also be improved by using more adjustable 
parameters; B, C, and ¢ can be adjusted independ- 
ently in all three configurations, thus increasing the 
number of adjustable parameters from 9 to 15. 
Theoretically, this should be done, as the effect of 
differences in these parameters is likely to have the 
same order of magnitude as the L(1+-1) correction.** 
The difficulty of carrying out least-squares calcula- 
tions with this number of parameters is great, even 
with SEAC. Some coding improvements and more 
experience in the behavior of the codes with this 
tvpe of problem are required. 
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Reflection and Transmission of Gamma Radiation by 


Barriers: 
Density Method ‘’ 


Monte Carlo Calculation by a Collision- 


Martin J. Berger 


Che collision density of photons in an infinite Compton-scattering medium was calcu- 


lated by random sampling 


The intensity, spectral composition, and angular distribution 


were then determined for radiation reflected and transmitted by plane-parallel barriers by 


carrving out appropriate integrations Over 


the collision density, 
absorbing properties of the medium and the effect of boundaries. 
done on the National Bureau of Standards automatic computer 


taking into account the 
All numerical work was 
SEAC). Sample results 


ire pres¢ nted for 0.66-Mevy radiation incident on water barriers of various thicknesses 


1. Introduction 


The Problem. The results of an investigation of 
the following boundary problem in’ gamma-ray 
diffusion are presented: A beam of monoenergetic 
gamma rays Is incident at a given angle on a plane- 
parallel barrier that has a finite thickness in one 
dimension but is infinite in the other two dimensions. 
What are the intensities and spectral compositions 
of the reflected and transmitted beam? 

Method of Solution. Ay analytical attack on this 
problem by the solution of the relevant transport 
equation leads to great mathematical difficulty. 
So far such an approach has been successfully carried 
out for homogenous infinitely extended media [1], 
but boundary problems remain unsolved. In order 
Lo bypass these difficulties, explore the problem, and 
obtain at least an approximate solution, this investi- 
gation employed the Monte Carlo (random sampling) 
method. All numerical work was carried out on the 
automatic computer of the National Bureau of Stand- 
ards (SEAC),. 

The calculation was divided, in regard to both 
method and execution, into two rather distinct parts: 
1) a stochastic part, in which photon random walks 
are generated by random sampling in an infinite 
nonabsorbing scatterer, whereby a collision density 
is obtained, and (2) a calculation of the reflected 
and transmitted radiation obtained by summations 
over collision densities, taking into account boundary 
conditions and absorption characteristics of the 
barrier. There were a number of reasons for this 
division. For one thing, it reduced the required 
computer memory capacity. More importantly, 
the random walks have a “universal” character, 
i. e., they can be used for the calculation of different 
boundary problems involving different geometries 
as well as different scattering and absorbing media. 
Not only will the repeated use of the same set of 


This work was supported by the Office of Naval Research and the Atomic 
Energy Commission Reactor Division 

Some of the results of this paper were presented at a symposium on Monte 
Carlo methods at the University of Florida in March 1954 


Figures in brackets indicate the literature references at the end of this paper 


random walks often lead to computing economy, 
but it may also increase the accuracy of a calculation 
Suppose, for example, that we wish to determine the 
difference in transmission for two beams incident on 
a barrier at different angles but otherwise identical. 
It will then be to our advantage to base the compari- 
son on the same set of random walks because 
irrelevant differences resulting from statistical fluc- 
tuations will thereby be largely eliminated. 


2. Gamma-Ray Random Walks in an In- 
finite Compton Scatterer 


2.1. Definitions 


The state S of a photon can be specified by a set 
of six quantities 


S=(10,¢,2,y,2), 


where / is the photon energy, @ and ¢ are angular 
coordinates describing the direction of motion (in a 
spherical coordinate system with the z-axis as polar 
axis), and v, y, and = are Cartesian coordinates of 
position. Let the state of the photon immediately 
after the nth seattering event occurring in a given 
random walk be denoted by S,(n=1,2,---), and 
let So denote the state in which it was introduced 
into the seattering medium. A random walk is 
then described by a sequence 


‘S. & Goace, 
each term (except So) depending stochastically on 
its immediate predecessor only (Markov 
The length, L, of such a sequence would be infinite, 
but in the present work the random walks are ter- 
minated when the energy, /, drops below 30 kev 
This arbitrary cutoff results in an average value 


pre CESS 


L)=i8 for an initial energy Ay=660 kev. A 
‘ The cutoff is justified on physical grounds because radiation at energies below 
1) kev is always so heavily absorbed that it makes only a negligible contribution 


to the emergent radiation flux 
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dom walks \Mionte Carlo 
he collision density resulting from the 


constitutes a 


clifi if radiation in an infinite homogenous 
( ton scatterel 
observed that for the problem of trans- 
na reflection by plane parallel barriers 
riables /, @ and one space variable, say 
ed Nevertheless it is worthwhile Lo 
the other three variables also The ran 
\ > are thus made applicable to other 


is with adiths rent geometrtl al cond! 


\loreove we can use. thre 


V prob 


sre set ol 


randol valks for barrier problems with different 
! ~ erence We alwavs set ‘iu QO and rotate 
he boundaries as required bv using, instead of cos 

~ os F COS a sit) ~ 

sit) y cos 
s the angle of ineider 
2.2. Random Sampling 

Next the sampling scheme ts described, by whic 
SHICCOSSIVE stutes S nre selec tea thre initial State 
being specified as S ..0.0.0.0.0 Calculations of 
st | ndom walks with this initial condition wil 
Serve the solution of problen s involving mono 
ner beams of energy /y and arbitrary direction 
| Ihe dence a Because ol the linearity ol tha 


possibli 1a) 


diff 
tutions for 


cumma-ray iIsiohn equation it is 


iacident 
energetl and angular composition hy superposition 


obtain al beams of speciliied 


of the results obtained with different /o’s and a’s 
Prion to the discussion of the detailed steps im thre 

sampling process, some comment ts in order concern 

ing the required random numbers. With a high 


speed computer the use of tables of random numbers 


is clearly mpractical lpstead, so-called pseudo- 


random numbers r,,(0<r,,<1) were used, and which 
were venerated ius required ith the eCourse of the 


Ceol }) Itution DY 


Dodd 


1! ethod cle Ve loped by () 
def ned by 


Taussky- 

Thev are the relations 

"RP . 
> 


R a R modulo 2 . R l. 


It mav be shown that the period of +, is 2°. 1. 
that a of 2 different numbers will be 
obtained repetition oceurs. Extensive test- 
ing carried out at the National Bureau of Standards 


Cape nee 


bye fore 


has shown that these pseudorandom numbers 
sutisiy the various accepted statistical criteria of 
randomness It Is an advantageous feature of this 


method that identical random walks 
created repeatedly for checking purposes, 
the initial number the 
recorded 

The 


vive'l 


inay be re- 
proy ided 


random used for walk is 


various the caleulation of S 


S. follow 


steps in 
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a. Energy Change 


The energy /,,, after a seattering is determined by 


7 - 
kh, A \dk kk), A \dE \==r, 3) 


Jt ‘ 


k(t.) is the 


per unit: path 


random number, and 
Klein-Nishina differential coetheient 
length) for Compton scattering with energy change 
from ki, to FB 

It was found convenient to carry out the random 
sampling for the wavelength \=me?//-, where me? is 


where / Is ii 


the electron rest mass, rather than for the energy 
directly If we make this change of variable. sub- 
stitute the explicit formula for the Klein-Nishina 
differential coefficient, and carry out the indicated 
Integrations, eq (5) Is transformed into 
( :, A \ \ 
l 2r,(A 1} log +r (X » 
/ \ 
\ N- ) 4 \ 2 
AX \ 1—2X,(A 1)| lo 
2x \/ \ 
: | Ka 2 1) 
2+A,)* } 
To use this equation to find X viven A, and +r, 
would be intolerably cumbersome, even when an 
automatic computer ts used Latter and Kahn [3 


numerically and presented a table 

for a grid of values in the range 
The tabulation has 
6 and used to obtain a polynomial 


evaluated Cg } 
of X nN 
Q<r<l, O.O5 <A, <10 
extended to X 
representation 


been 


] J 
\ . >. Add.) 5 


This required extensive bivariate interpolation, 
which was carried out by a convenient matrix tech- 
nique. This technique is described in the appendix 
because it does not appear in the standard literature 
on numerical analysis 

The following coefficient-matrices, <1,,, were found 
to give representations of A accurate to 1 percent 


" 
s 44 . sur 
6u x2 
! Jn 4 4 
j 27s ~ 
\ 
i 
shin os s. 4 ; 
774 1. ol 1 SH s 1. 330 
au THO 01 "4 0. 1780 
“ ON 0) ” ; ’ 
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b. Change of Direction 


Deofl ection and wavelength shift are related by the 
Compton equation 

(‘ts WwW | Xr 

where w, is the angle between the directions of motton 

immediately before and after the n+ 1st seattering 

The dire tion of the scattered photon depends on 

w, and on an azimuth angele x,, which is a random 


quantity distributed uniformly between 0 and 27 
We need the sine and cosine of x, rather than x 
itself Hence we take advantage of the following 


convenient computational scheme suggested by Von 
Neumann {4 random numbers a and b 
satisfying the condition a?+6?<1, and let ¢ be a 
random number that is equal to +1 with probability 


( ‘hoose 





12. Then 
ahe ad h- - 
COs X and sin x : ‘ 
ad h ad h 
Krom the deflections w, and x,, the new angular 
coordinates are then determined by the following 
trigonometric relationships: 
cos 4 cos 7, COS @ t+-siné Sin w, COS X 7 
Sin xX sinw 
Sill (@ oO 
sin 6 > S 
COS @ cos 7, cos 6 
COS (Od @ : : 
sin @. sin @ J 
c. Displacement 


The position of the (n+ 1)st seattering is defined by 





. 
sin 4, cos ¢ 
/ / - oY 7] 
sin @, sin ¢, 
7 logr > ) 
Li (F,) 
cos @ 
. Og 7 
uu E,) J 


where ris a random number, and 


*f 


K(k, Ajdk 


3. The Transmission-Reflection Boundary 


Problem 


3.1. Summation Over the Collision Density 


The information obtained by sampling J photon 
random walks is contained in the set of state vectors 
S,(p(n=0,1,- ++ Ly; 7=1,2,- ++ J), which describe the 
state of the jth photon immediately after its nth 
scattering. We consider form a 
of representative points (a collision 
density) mn a six-dimensional collision-density 
The solution of boundary problems can be obtained 
by performing appropriate sums over the collision 
density. 

In the reflection-transmission problem, one con- 
siders a plane-parallel barrier between the planes 

Oand z=¢t. A beam of photons is assumed to be 
incident on the face 0) 

Reflected and transmitted radiation will be classi- 
fied according to energy and direction.’ Let R 
denote the average number of photons per incident 
photon that are reflected from the face QO with 
energies in the 7th interval and directions in the Ath 
interval. Similarly, let 7,, denote the average num- 
ber of photons transmitted through the face f 


vectors to 


sample 


these 
set 


space 


R, and Ty are obtained by summations over the 
sample collision density as follows 
IA , , ' 
Ri, pou 2B pP ()Q,9,O) *()). 10) 
The factors in this formula have the following 


meaning: 

(a) The factor B,()) takes into account the bound- 
ery conditions. It is equal to unity when the in- 
equalities 0O<2,,(7) <¢ are satisfied for m=0,1, - ++ n, 
and zero otherwise. In other words, a contribution 
to the reflected radiation is obtained only if the 
photon in the under consideration has not 


previously gone out of the barric 
-n(D) } 11) 


(b) The factor 
P,(j)= I exp {- SeC Om\J) I, 

= UT, walEn (3) 
is the probability that the photon has not 
absorbed in the jth history prior to reaching state 
SQ u, Is the absorption coefficient. ) 

(c) The quantity {,(j,2°) denotes the probability 
that a photon starting from state S,(7) will cross the 
plane ‘ without further scattering or absorption. 
It is defined as follows: 


state 


1(J) 


been 


a; vanishes except when sec 6,(j)){2’ nn) 

>0, in which case 
° , sec 6,,()) . -_ ; 
@,(j,7) exp <4 — (K(; [2 =n())I 12) 
ul £,())| 

Che method described in this section is closely related to suggest 
H. Kahn, [3 

‘2 energy intervals and 10 angular inter have been used 
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s the total linear attenuation coefficient 
is the energy-angle classification factor 


to unity when /,(7) is in the 7th energy 

and @ in the Ath angular interval, and 
therwise 

msmitted radiation 7’, is calculated by a fol 


i similar to eq 10) but with ‘J, 7,0 replaced by 
We shall use the collision density also for the eval- 
of the flux in an infinite medium in which a 
source is embedded in the plane 0. By 
comparing the results for this problem with those 
for the barrier, the nature of the boundary effect can 
he established \loreover, the Monte Carlo results 
for the infinite medium can be checked against those 
obteimed by systematic moment-method calculations 
Such a comparison is of use in establishing the valid- 
it\ and accuracy of the random-sampling approach 
Let By(2’) denote the average number of photons 
crossing the plane from both directions (our 
main interest will be in the planes that form the 
barrier boundaries Oand 2’=—/ The necessary 
over the collision density is similar to 
except. that the boundary factor 


ili}? 


rizbiie 


summation 

formula (10 

B j) is omitted 
When the sampled collision density is used for the 


simultaneous determination of P?,,, 7... B,.(0), and 
B,(t), all states S,(7) make a contribution to the 
final score except those for which ) 0 and 
cos 6,(7)< OQ. or } tf and cos é J 0) This is in 
contrast to a direct stochastic analogue method, 


which would allow only those sections of the random 
walks containing an actual crossing of the boundaries 
to contribute to the The more 
gives proper eredit for “near 
i. € collisions taking place close to a 
boundary, for which the factor ‘y is very 
to unity 


score elaborate 


scoring procedure 
misses,” 


close 


3.2. Fluxes and Buildup Factors 


The SEAC calculation is designed to vield the 
quantities Ry, Ty, and B,,, which represent the 
average number of photons crossing the planes 0) 
and t per unit area (the number flux Radi- 
ation detectors often do not measure the number 
flux directly, so that it is desirable to calculate 
various related types of flux These can all be 
obtained from the number flux through multi- 
plication by suitable weight factors, depending on 
the energy and direction of the photons kor 


example, multiplication of the number flux by the 
average energy, (/:,), vields the energy flux, division 
by the value, value 
vields fluxes through a unit area whose normal is in 
the direction of motion of the radiation This is 
the flux that a nondirectional detector would 
measure. Some inaccuracy is introduced into these 
flux COnVEerSLONS the averages Ie and 
extent arbitrary. Arithmetic 


average absolute cos 6 


because 
cos are to some 


means have been used 


According to the common in. shielding 
calculations the results for the total flux (integrated 
over-all directions and spectral energies ) are pre- 
sented in the form of buildup factors, defined as 
the ratio of the total to the unscattered flux [5]. 
For the various types of flux described above, the 
buildup factors for an infinite homogenous medium 
have the following form: 


usage 


Number buildup factors 


B, > Ba 


Pe eX) | My = COS a@ 
. B eXp Unt COS Qa 
By=> | 13 
— cos # COS a 
Energy buildup factors 
Bb, >> / Bf CXp My. COS Qa 
: B 
i, > [ exp Unt COS @), 14 
cos # COS a 
where wo is the total attenuation coefficient of the 
source radiation, and ¢ is the distance between the 


source plane and the plane of observation. — In 
order to obtain the corresponding buildup factors 
Ty. Ty, Te, and Ty, one replaces By in eq (13) 
and (14) by T and by the barrier thickness, ¢. 
The reflection buildup factors lew. Ry. Rx. and R: 


are obtained by using /??,, and setting 7=0 


4. Results’ 


4.1. Energy Spectra 


In figures 1, a, b, ¢, and d, energy spectra are 
shown in histogram form for the number flux of 
photons reflected and transmitted by barriers, 
reflected by a semi-infinite medium, and for the 


flux in an infinite medium through planes that 
correspond to the barrier boundary planes. The 
histograms pertain to scattered radiation only. 
Spectra are shown for normal incidences (a=0 


The oblique bar- 
») 


and for oblique incidence a—b60 
rier thickness is the same for both cases: y 
The shaded areas between the “finite slab” and the 
“infinite medium” histograms are a measure of the 
boundary effect. It can be seen that the boundaries 
greatly reduce the amount of soft radiation below 
200 kev. The reason for this is that in an infinite 
medium soft radiation is mainly due to photons that 
have overshot the boundary plane and have been 
turned around in a collision in the close vicinity of 
the boundary, making their second passage through 
the boundary with resultant low energy. 


uf COS a 


Based on the analy of 400 random walks for each problen 
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4.2. Angular Distributions 


Figure 2, a. shows the angular distribution of the 
reflected number flux from a barrier with oblique 
thickness pot =1, and from a semi-infinite medium 
for obliquity angles of incidence a=0° and a=60 
Figure 2, b, shows the angular spectra of transmitted 
scattered photons for normal incidence and barrier 
thickness pof=1 and pof =—4. 


mn the barrier (or rele 


ised within a semi-infinite medium) with obliquity @ 
6°; ¢, Transmission: obliquity a=0° (normal incidence); d, Transmission; obliquit 


4.3. Buildup Factors 


The various buildup factors for transmission, 
reflection, and an infinite medium (as defined in 
section 3.1 and 3.2) are listed in table 1 for the 
number flux, and in table 2 for the energy flux for 
(.66-Mev radiation in water. Estimated standard 
deviations are indicated, which were obtained by 
dividing the 400 photon histories into groups of 50 
each and computing the dispersion of the buildup 
factors obtained for the various groups. 
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5. Discussion 


The results presented have a preliminary charac- 
ter. Greater sample sizes are desirable, as well as 
the extension of the calculation to other energies 
and media. Moreover, the scoring method used its 
not expected to be very successful for penetrations 
of depth of >4 because excessively large sample 
sizes would be required. Prior to extending the 
calculations by the collision-density method, other 
Monte Carlo approaches have been explored, in- 
cluding a semianalytic method designed for dealing 
with deep as well as shallow penetrations [6]. Table 
3 gives a comparison of the buildup factors obtained 
in this investigation (a) with those found by the semi- 
analytical Monte Carlo method, (b) with the Monte 


Carlo results of a group at the Naval Research 
Laboratory |7], and (c¢), for an infinite medium, with 
calculations carried out according to the moment 


method of Spencer and Fano |1, 8] The agreement 
is generally indicating that the shallow 
penetration of gamma radiation, in the presence of 


vood. 


boundaries, can be calculated accurately by the 
| poe 38 (‘ompa son of bu dup factors obtained in arious 
( P 0.66-Mev radiation 7) 
t nalvt NiIont Car 1 t} r : Nlonte 
KI i ‘ I ent method [s 
4 
| 
xu 
21 
4 
1.2 
2. 4s } 2.4 
07 ; Z su 
~ ri 7 
1 or 
2. | 2 0 
s , 2 
21 7.82 
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| Do p: : p Pi Po 
Pit Po) 


Pr pi 
Po Px) (Po 
L Pr Ps 
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Po) 
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Pi) (Po 





P2) 


Monte Carlo method. The results for Ay and RF, 


at 0.66 Mev are consistent with the results of a 
Monte Carlo calculation at 1 Mev by Hayward and 


Hubbell [9] 
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7. Appendix. A Method of Bivariate Inter- 


polation 
Suppose we know a bivariate function f(p,q) at 
Mf+1)(N+1) gridpoints (ps,g:) (s=0,1, -, M; 
t=0,1, ...,N). The problem is to use this knowl- 


edge for determining the expansion coefficients in a 


polynomial expansion 


V/ \ 
M(p, M= Do DS Avp'y’. 7) 
Let F=f(p..q.) be the known matrix of functional 
values, and A= (A the desired coefficient matrix. 
From the variables p and gq we form the matrices 
P=(p.)' and Q=(q,)'. In matrix form, eq (17 
becomes 

F=PAQ’, i 

where Q’ is the transpose of Y@. Hence 
A=-P"'Pig" ty’. 19 


P and @ are so-called alternant matrices whose in- 
verse is well known [10]. The prescription for ob- 
taining the inverse is simplest to describe by means of 
anexample. When / has rank 3 
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(Po P2) (Pot P1) on 
(Pi-— Po) (Pi— P2) = (Pe Po) (Po— Pr) os 
l l 
(Pi— Po) (Pi— P2) = (Po— Po) (P2 ri)J 
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The sum of the elements in the first row of each 
inverse matrix is 1, and the sum of the elements in 
each other row is zero. This is also true for the case 
of nonequidistant grid points, and provides a useful 
arithmetical check 
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